o 


Cruft  Laboratory 


CJ>  o 

CD  —I 

CD  ^ 


o 


CQ 

Q 

LlJ 

C.D 

O 


CD 


<c 

CD 


CO 

<c: 


Harvard  University  »  Cambridge,  Massachusetts 


THEORETICAL  AND  EXPERIMENTAL  STUDIES 
OF  ANTENNAS  AND  ARRAYS  IN 
A  PARALLEL -PLATE  REGION 


October  5,1962 


Scientific  Report  No.  15  (Series  2) 

Contract  AF  19  (604)  -  4118  Project  5635,  Task  563502 
AFCRL- 63-90 


I  -  - -  1 

Prepared  for 

Electronics  Research  Directorate 
Air  Force  Cambridge  Research  Laboratories 
Office  of  Aerospace  Research 
United  States  Air  Force  •  Bedford,  Massachusetts 


AFCRL-63-90 


THEORETICAL  AND  EXPERIMENTAL  STUDIES  OF  ANTENNAS 
AND  ARRAYS  IN  A  PARALLEL-PLATE  REGION 


by 


B  .  Rama  Rao 


Division  of  Engineering  and  Applied  Physics 
Cruft  Laboratory 
'  Harvard  University 
Cambridge,  Massachusetts 


Contract  AF19(604)-4118 
Project  5635  -  Task  563502 
Scientific  Report  No.  15  (Series  2) 
October  5,  1962 


Prepared  for 

Electronics  Research  Directorate 
Air  Force  Cambridge  Research  Laboratories 
Office  of  Aerospace  Research 
United  States  Air  Force 
Bedford,  Massachusetts 


Requests  for  additional  copies  by  Agencies  of  the  department  of  Defense, 

their  contractors,  and  other  Government  agencies  should  be  directed  to  the: 

ARMED  SERVICES  TECHNICAL  INFORMATION  AGENCY 
ARLINGTON  HALL  STATION 
ARLINGTON  12,  VIRGINIA 

Department  of  Defense  contractors  must  be  established  for  ASTIA  services,  or 
have  their  "need-to-know"  certified  by  the  cognizant  military  agency  of  their 
project  or  contract. 


All  other  persons  and  organizations  should  apply  to  the: 

U.  S.  DEPARTMENT  OF  COMMERCE 
OFFICE  OF  TECHNICAL  SERVICES 
WASHINGTON  25,  D.  C. 


SRI  5 


.  ** 


THEORETICAL  AND  EXPERIMENTAL  STUDIES  OF  ANTENNAS 
AND  ARRAYS  IN  A  PARALLEL-PLATE  REGION 

PART  I  -  THE  CURRENT  DISTRIBUTION  AND  IMPEDANCE  OF  AN  ANTENNA 
IN  A  PARALLEL- PLATE  REGION 

by 

•  B,  Rama  Rao 

ABSTRACT 

The  integral  equation  for  the  current  on  the  antenna  has  been 
formulated  using  Green's  theorem  and  a  suitable  Green's  function  based 
on  the  method  of  images;  the  resulting  equation  has  been  solved  by  expressing 
the  current  as  a  Fourier  series  in  terms  of  waveguide  modes.  It  has  been 
shown  that  in  the  case  of  thick  antennas  the  logarithmic  term  in  the  current, 
due  to  the  idealized  delta-function  generator,  makes  a  marked  contribution 
over  a  significant  length  of  the  antenna  and  the  so-called  'transition'  region 
is  much  larger  than  for  the  free-space  antenna.  The  admittance  of  the 
theoretical  model  can  be  brought  into  good  agreement  with  the  experimental 
results  by  subtracting  out  the  infinite  gap  capacitance  and  replacing  it  by 
an  empirical  correction  term  to  account  for  the  input  susceptance  of  the 
actual  feeding  gap.  An  anomaly  in  the  behavior  of  the  antenna  at 
resonance,  noticed  by  Lewin,  can  be  avoided  by  assuming  a  complex 
propagation  constant  in  the  guide,  due  to  small  but  finite  losses  in  the  image 
planes.  Closed  form  expressions  for  the  current  and  admittance  have 
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been  obtained  and  compared  with  experimental  results.  The  general 
agreement  is  found  to  be  good. 

To  avoid  complications  due  to  the  feed-point  singularity,  the  admit¬ 
tance  has  also  been  obtained  by  a  variational  method.  There  is  good 
agreement  with  experimental  results  for  antenna  lengths  in  the  range 
0  <  kQh  <  3.  7  ;  but  the  theory  deteriorates  for  higher  values  of  kQh. 

Extensive  experimental  measurements  have  been  made  on  the  antenna 
and  a  description  of  the  apparatus  is  included  in  the  report. 


SRI  5 


-3- 


1.  INTRODUCTION 

Much  attention  has  been  focused  during  recent  years  on  obtaining 
mathematically  rigorous  solutions  to  a  large  variety  of  problems  concerning 
antennas  radiating  in  free  space  (1);  depending  on  the  complexity  and  nature 
of  the  problem  the  techniques  employed  range  from  rigorous  Fourier- 
transform  and  Fourier-series  solutions  to  the  more  approximate  iterative 
.and  variational  techniques.  In  most  cases  the  theoretical  results  have 
been  supported  by  comprehensive  experimental  investigations.  In  contrast, 
however,  the  general  problem  of  determining  the  current  distribution  and 
driving-point  impedance  of  antennas  radiating  into  waveguides  has  received 
scant  attention  until  quite  recently.  Most  of  the  early  investigations  (2) 
were  confined  to  obtaining  the  impedance  of  small,  'filamentary'  antennas 
employing  the  emf  method  with  an  assumed  current  distribution  on  the 
antenna.  The  inherent  inadequacy  of  such  a  method  as  compared  to  the  more 
sophisticated  integral  equation  method  and  its  variational  modification  has 
been  pointed  out  by  King  (1).  In  practice, it  is  also  not  permissible  to  lay 
such  severe  restrictions  on  the  finite  thickness  and  size  of  the  antenna. 
Lew'in(3>4)  has  lately  given  the  problem  a  more  careful  mathematical 
treatment.  In  a  recent  paper  (3),  he  has  considered  the  radiation  from  a 
linear  antenna  in  a  parallel-plate  region  using  a  modal  expansion  method. 

The  solutions  he  obtains  for  the  current  and  input  admittance  are  not  in  a 
closed  form  and  cannot  be  evaluated  numerically  for  comparison.with 
experimental  results.  Furthermore,  he  has  noticed  certain  anomalies 
concerning  the  behavior  of  the  antenna  near  resonance  .which  cannot  be 
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explained  physically;  all  of  these  aspects  deserve  further  examination,  both 
theoretically  and  experimentally.  The  main  motivation  for  the  present 
investigation,  however,  arose  while  considering  the  propagation  of  surface 
waves  along  an  antenna  array  placed  in  a  parallel-plate  region;  for  treating 
this  problem  rigorously  it  was  necessary  to  know  the  current  distribution 
on  the  individual  elements  of  the  array  and  the  mutual  coupling  between  the 
elements.  An  exact  solution  to  the  simple  problem  of  a  single  antenna  in 
the  parallel -plate  region  is  a  necessary  prerequisite  towards  obtaining  a 
better  understanding  of  this  more  complicated  problem. 

In  this  report  an  integral  equation  technique  has  been  employed  to 
obtain  the  current  distribution  and  admittance  of  the  antenna.  The  mathematical 
procedure  adopted  here  closely  follows  that  used  for  solving  the  free-space 
antenna  problem.  Since  these  two  cases  are  twin  aspects  of  the  general 
problem  of  antenna  radiation,  a  comparison  between  them  might  prove  to  be 
instructive.  Another  attractive  advantage  of  the  integral  equation  method  is 
that  it  can  readily  be  extended,  as  outlined  by  King  (1)^  for  obtaining  the  current 
distribution  on  coupled  antennas  or  in  an  array  of  antennas, 

To  facilitate  experimental  measurements  the  antenna  was  driven  near 
its  base  by  means  of  a  coaxial  line  connected  to  one  of  the  ground  planes; 
this  is  in  contrast  to  Lewin1  s  case  (3)  where  the  antenna  was  driven  at  its 
center  between  two  plates,  but  it  was  not  a  practically  feasible  scheme  and 
had  to  be  discarded.  Other  than  this  slight  difference,  the  problem  treated 
here  and  the  one  considered  by  Lewin  are  identical. 
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In  Sec.  2  the  integral  equation  for  the  current  has  been  solved  by- 
expressing  the  current  on  the  antenna  as  a  Fourier  series  in  terms  of  wave¬ 
guide  modes.  Special  attention  has  been  paid  to  the  following  aspects  of 
the  problem:  (1)  the  logarithmic  singularity  in  the  current  and  its  effect  on 
admittance,  (2)  the  behavior  of  the  antenna  at  resonance,  (3)  current  on 
an  infinitesimally  thin  antenna,  (4)  the  charge  distribution  on  the  antenna, 

(5)  the  effect  of  assuming  a  filamentary  kernel  in  the  integral  equation.  Closed- 
form  solutions  for  the  current  and  admittance  have  been  obtained  and  compared 
with  the  experimental  results.  The  general  agreement  is  very  good. 

While  the  Fourier-series  method  provides  accurate  solutions  for  the 
current  distribution  and  conductance  of  the  antenna,  the  susceptance  values 
obtained  from  it  are  rather  ambiguous  owing  to  the  infinite  gap  capacitance  of 
the  delta -function  generator.  In  Sec.  3  the  impedance  of  the  antenna  has 
been  obtained  by  a  variational  technique  using  a  suitable  trial  function  for  the 
current.  The  variational  method  ignores  the  singularity  at  the  driving  point 
and  provides  an  approximate  but  continuous  solution  for  the  current  distribution 
over  the  entire  length  of  the  antenna.  The  theoretical  and  experimental 
admittance  values  are  in  excellent  agreement  until  kQh  =  3.7  but  deteriorate 
rapidly  thereafter. 

Section  4  is  devoted  to  experimental  phase  of  the  research.  A 
description  .is  given  of  the  apparatus  and  experimental  techniques  used  for 
determining  the  current  distribution  and  admittance  of  the  antenna. 
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2.  FOURIER-SERIES  SOLUTION 
Formulation  of  the  Problem 

The  antenna  under  consideration  is  the  center  conductor  of  the  coaxial 
cable  which  extends  inside  the  parallel -plate  region  and  has  its  upper  end 
connected  to  the  top  plate.  Fig.  1.- 1  shows  the  geometry  of  the  idealized  model 
driven  by  a  delta -function  generator.  A  cylindrical  coordinate  system 
(r  ,  0  ,  z)  is  introduced  and  the  axis  of  the  antenna  of  radius  a  is  made  to 
coincide  with  the  z  axis.  The  two  infinite,  perfectly  conducting  image  planes 
are  situated  at  z  =  0  and  z  =  h.  The  antenna  is  assumed  to  be  excited  by 
a  slice  generator  situated  at  z  =  _+  5  .  The  equivalent  structure  obtained 
by  considering  successive  reflections  in  the  two  planes  is  that  of  an  infinitely 
long  antenna  fed  by  generators  along  the  entire  length  spaced  at  intervals 
depending  on  the  distance  h  between  the  two  planes.  Alternatively, the 
antenna  may  also  be  considered  as  an  infinitely  long  collinear  antenna  array 
of  half  length  h  ,  placed  end  to  end. 

Because  of  angular  symmetry,  the  only  nonvanishing  components  of  the 
field  are  E^,  Er  and  Hg  .  Since  E^  vanishes  at  the  two  plates,  the  z  component 
of  the  vector  potential  Az  satisfies  the  Neumann  boundary  condition 


8A_ 


8  z 


=  0  at  z  =  0  and  h 


(l.D 


The  general  equation  satisfied  by  A  at  all  points  outside  the  conductors  is 


(V2  +  k  2)  A  = 
o  z 


(1.  2) 


where  kQ  =  2  IT  /X.  is  the  wave  number. 
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FIG.  1-2.  Image  Representation  of  Greens  Function 


FIG.  1-3  Application  of  Greens  Symmetric 
Theorem  to  Antenna  Configuration. 
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The. solution  of  this  equation  may  be  obtained  by  applying  Green's 
symmetric  theorem 


J*  j^uV2  v  -  yV2  u  dr1  =  J* 


UJI _ 8u 


8n 


8n 


d<r' 


(1.3) 


with  u  =  and  v  =  G,  where  G  is  the  Green's  function  satisfying  an 
inhomogeneous  wave  equation. with  periodically  distributed  sources  (see  ensuing 
discussion) 


V2 


+  k 


G  = 


I 


m 


5  (r  -  r  ') 
m 


(1.4) 


and  the  following  boundary  conditions: 

(1)  G  is  continuous  except  at  r  =  r  '  where  it  has  a  logarithmic  singularity, 

(2)  G  satisfies  the  radiation  condition  at  r  •»  oo  and  the  Neumann  boundary 


condition 


8G 


8  z 


=  0  at  z  =  0  and  h. 


k  is  the  free-space.wave  number:  r  and  r  1  are  the  vector 
o  r  ’  m 

coordinates  of  the  field  and  source  points  respectively;  6  (r  -  rm')  is  the  three- 
dimensional  Dirac  delta  function  defined  in  the  usual  manner  as 


S  (  r  -  r  ')  =  0 


f 


5  (r 


m 

')  dr  =  1 


if  r  includes  r  1 
m 


=  0 


if  r  excludes  r  1 
m 


Image  Representation  of  the  Green's  Function 

The  appropriate  Green's  function  for  the  problem  can  be  constructed  by 
using  the  principle  of  images  (see  Fig.  1.  2).  Let  the  primary  source  point 


SRI  5 


-7- 


The  solution  of  this  equation  may  be  obtained  by  applying  Green's 
symmetric  theorem 


r 

uV2  v  -  vV2  u 

dr'  =  f 

C  SV 

•  -v  8u 

d 

8n 

8n 

T 

-  - 

Z 

_ 

da' 


(1.  3) 


with  u  =  and  v  =  G,  where  G  is  the  Green's  function  satisfying  an 
inhomogeneous  wave  equation  with  periodically  distributed  sources  (see  ensuing 
discussion) 


v2 


+  k 


G  = 


6  (r  -  r  ') 
m 


(1.4) 


and  the  following  boundary  conditions: 

(1)  G  is  continuous  except  at  r  =  t^'  where  it  has  a  logarithmic  singularity, 

(2)  G  satisfies  the  radiation  condition  at  r  -*■  oo  and  the  Neumann  boundary 


condition 


8G 


8  z 


0  at  z  =  0  and  h. 


k  is  the  free-space  wave  number;  r  and  r  '  are  the  vector 
o  r  ’  m 

coordinates  of  the  field  and  source  points  respectively;  6  (r  -  r  ')  is  the  three- 
dimensional  Dirac  delta  function  defined  in  the  usual  manner  as 


6  (r  -  r  ')  =  0  r  / 


f 


6  (r  -  r  ')  dT  =  1  if  t  includes  r  1 

m  m 


=  0 


if  t  excludes  r  ' 
m 


Image  Representation  of  the  Green's  Function 

The  appropriate  Green's  function  for  the  problem  can  be  constructed  by 
using  the  principle  of  images  (see  Fig.  1.  2).  Let  the  primary  source  point 
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be  situated  at  z'.  With  consideration  of  its  successive  reflections  in. the  two 
planes  an  infinite  series  of  image  sources  is  obtained.  These  are  located 
at  2mh-z'  where  m  =  0,  +1,  +  2.  .  .  and  at  2m h  +  z1  where  m  =  +  1,  +  2, . . 
as  shown  in  Fig.  1.  2.  The  resultant  effect  of  all  these  image  sources 
is  such  that  G  satisfies  the  Neumann  boundary  condition  8G/8z,  =  0 
at  z  =  0  and  h. 

In  view  of  the  above  discussion  and  in  analogy  with  free-space  antenna 
theory,  a  suitable  Green's  function  satisfying  all  the  boundary  conditions 
mentioned  in  Eq.  1.  4  is 


-jk  R 
J  o  m 


linr 


(1.5) 


where  Rm  =  |r  -  r _ '  I  denotes  the  vector  distance  between  the  field  and 

m  1  m 

source  points. 

Now  consider  Eqs,  1,  2  and  1.  4  and  apply  Green's  symmetric  theorem 
to  a  pill-box  shaped  surface  of  the  form  shown  in  Fig.  1.  3.  It  follows  from 


Eq.  1.  3  that 


s 


G  (V2  A  +  k2A)-A  (^G  +  k  2  G) 


z  o  z 


dT* 


S1*WS4 


|G  «A  /bn  -  A  — 

\  z'  z  Bn 


do-'  (1.6) 


where  Sj  is  the  outer  surface  of  the  pill-box,  and  Sj  are  the  top  and  bottom 
surfaces  and  the  circumference  of  the  antenna.  As  recedes  to  infinity, 
the  integral  over.S^  reduces  to  zero  since  . both  A  and  G  vanish  at  infinity  due 
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to  the  radiation  condition.  Also  the  integrals  over  and  Sj  contribute 

9  G 

nothing  since  both  8  A^/an  and  —  vanish  on  the  image  plates. 

Furthermore,  since  A^  =  Ag  =  8A/86  =  0,  and  since  the  method 

of  excitation  is  such  that  only  the  z  component  of  the  current  exists 
on  the  antenna,  has  to  satisfy  the  following  boundary  condition  on  the 
surface  S  .  of  the  antenna. 


where  is  the  current  density  on  the  surface  of  the  antenna. 

Thus,  of  the  four  functions  in  equation  1.  6  only  dA^/dz  is 
discontinuous  across  S^.  Hence, 


K  (z * )  Gdo-' 
z  '  ' 


(1.8) 


where  r  1  refers  to  the  vector  coordinates  of  the  primary  source  point. 
With  I(z')  =  2ffa  Kz(z*),  Eq.  1.8  becomes 


27 T 


A  (r)  =  — — 
z  v 


I  (z1)  dz' 
z  ' 


dj( 


2tt 


(1.  9) 


o  o 

where  /  =  61  -  0  and  G  is  given  by  1.  5. 

From  antenna  theory  the  distance  Rm  in  the  kernel  is  given  by 

1/2 


R 


(  z  +  z1  -  2mh)^  +  (2a  sin  0’/2)  ^ 


(1.10) 
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where  K  is  the  modified  Bessel  function  of  the  second  kind.  With  1.14  and 
o 

1.12  in  1.11  it  follows  that 


■jko 

e 

(z-z1  -  2mh)^  +  (2  a  sin  j/  /2)^ 

' _  _ i 

1 1/2 

(z-z1  -  2mh)2  +  (2asin  j //2)2 

1/2 

OD 

■I  ± 

m=-oo 

Similarly, 


j  mir  (z-z1) 

IT 


K  J(2asin-4-)2  lJ%7~)2 -K‘ 2 


(1.15) 


_jko  ^ 


z  +  z'  -  2mh)2  +  (2  a  sin  <jt  / 2 )2 


(z  +  z1  -  2mh)2  +  (2a  sin  /2)2 


j  mrr  (z  +  z') 
h 


K  /  (2a  sin  ) 


-  kQ  =  rm  is  easily  recognized  to 


-1  T  «  K  Juaain  4-)2  (^l)2-ko2  .  ( 

m=  -oo  _ L _  J 

In  the  above  equations  J  ~  kQ2  =  is  easily  recognized 

be  the  mode  propagation  constant  in  the  parallel -plate  guide. 

It  is  now  possible  to  add  (1. 15)  and  (1. 16)  and  combine  the  terms 
symmetrically  with  respect  to  m  =  0  .  With  this  result  in  1.  11  and  the 
substitution  of  2  0  for  $  ,  the  following  result  is  obtained 

h  oo  f  77  'j 

Az(r)  =  -25TT  f  I(z')dz'  "5rTT  I  Ko(2koa  V8in  e  )  d0f 

o  m=-oo  Lo  J 


(1.16) 


cos  mirz  cos  mir  z' 

"i r  k 


(1. 17) 
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(E  e)  =  -  V  6  (z) 
z  r=a  ° 

Hence,  the  electric  field  at  any  point  along  the  antenna  is 

(E  )  =  Z'lfz)  -  V  6  (z)  0  <  z  <  h 

z  o 

r=a 

For  a  perfectly  conducting  antenna  Z1  =  0  and 

(E  )  =  -  V  6  (z) 

z  o 


(1.  22) 


From  1.22  and  '1.20 


Vq  6  (a)  =  jtQ  f  K{.z,  a.)  I  (a')  dz< 


(1.  23) 


where  the  kernel  of  the  integral  equation  K(z,  z')  is  given  by 

7 r 

/  Ko«2ko  a  sin  0  )  d  0 


K(zj  z1)  —  - ^ — 

2tT  h 


oo  7T 


-2  I  /fo<2ko^mSine4  d0 


m=  1  o 


l  do 

1 

; 

Ll  °) 

m7T  z  m77  z  ' 

cos — r —  cos  —r. -  . 

h  h 


(1.  24) 


The  integrals  appearing  in  1.  24  can  now  be  evaluated  by  means  of  the 
appropriate  integral  representation  for  the  Bessel  functions.  The  formulas 
used  in  the  derivation  below  are  all  obtained  from  Magnus  and  Oberhettinger  (5), 
which  will  hereafterwards  be  briefly  referred  to  as  (M  and  O). 
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Cons  ider  the  first  integral 


M  ‘  II 

P  =  J1  Ko(akQaYm  sin  0)  d0  =  KQ(Cjn  sin  0)  dO  (1.  25) 


where 


c  =  2k  av 
m  o  m 


The  integral  representation  for  the  Kq  function  is  given  by  M  and  O  (5) 
page  27  as 


Ko(ax)  = 


P  cos  X  t 

y  /  2  2 

o  ./a  +  t 


dt  (a,  x  real  and  positive  )  . 


With  this  result  it  follows  that 


IT  oo 


■<f  S 


[t  sin"o1 


dt  dO  = 


w  v  /  2  2 

o  o  /  c  +  t 

v/  m 


oo  IT  _ 

cos  t  sin  ol  < 

o  /c^  +  t^o  —  -> 

\J  m 


“  J  (t) 

=  IT  P  ■  — ■  ■  °— 

c/  / 


dt.  (1.26) 


V  /  2  .  .2 

o  /  c  +  t 

v  m 


Use  .is  now  made  of  Watson's  formula  (5)  [M  and  O  p.  30] 


yu 

!v/2  (  T  ax>  Kv/2  (  T  ax)  =  ,/ 


J  (at) 
v 


o  t“  + 


(Re  v  >  -  1,  a  real  and  positive,  |arg  x  |  <ir  /Z)  . 


With  v  =  0  in  the  above  equation  and  with  the  substitution  of  this  result  in 


1.  26  ,  it  follows  that 
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=  / 


K  (2k  a  y  sin  0  )  dO  =  7T  I 
o  o  'm  c 


(k  a.y  ) 
o  'm 


K 


(k  a.y  ) 
o  'm 


(1.  27) 


where  I  and  K  are  the  modified  Bessel  functions  of  the  first  and  second  kind, 
o  o 

For  imaginary  arguments  the  identities  Iq  (jz)  =  Jq  (z)  and 

/  V 

Kq  (jz)  =  -j  ~  Hq  (zj  apply, where  Jq  (z)  is  the  Bessel  function  of  the 

(2) 

first  kind  and  H  '  '  (z)  is  the  Hankel  function  of  the  second  kind.  With  m  =  0 
o 

in  1.  27  and  with  these  several  identities  it  follows  that 


f  Ko  «2ko 


a  sin  0)  dO  =  - 


jTT 


II 

f  Ho 


(2) 


(k  a  sin  0)  dO 
o 


2 


^ —  J„  (k„a)  H^(2)  (k„a)  .  (1.  28) 


o  o  o  o 

(See  Appendix  I  for  rigorous  derivation  of  this  formula). 

The  substitution  of  1.  27  and  1.  28  in  1.  24  leads  to 


-k 


K(z,  z1 


2nh 


00 


+  2  y  1  (k  )  K  (k  a.y  ) 
l_j  o  o  m  o  o  m 

m=l 


j  _I_  J  (k  a)  H  (2)  (k  a) 
J  2  o  o  o  '  o 

.2  “ 


mTT 

k  h 

,  o 


1 


cosm7Tz  cosmirz1 
h  Ki 


for  0  <  kQh  <  7T  .  (1.  29) 

When  the  separation  of  the  parallel  plate  exceeds  a  half-wavelength, 
that  is  when  7T /kQh  <  1,  the  term  corresponding  to  m  =  1  has  to  be  modified 
to  account  for  the  propagation  of  the  first  higher-order  mode  in  the  parallel- 
plate  medium.  The  corresponding  expression  for  the  kernel  is 
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The  current  on  the  antenna  may  now  be  expressed  in  terms  of  the 
Fourier  series 

oo 

I(z)  =  1+2  )  I  cos  m7T  /h  . 

m=l 


(1.  33) 


After  the  substitution  of  1.  33  in  1.  31  the  integration  may  be  performed  to  obtain 

I  =  V  /]%  h2  a 
m  o'Jso  m 

The  current  distribution  on  the  antenna  is  then  given  by 

V 

I(z)  =  I"(z)  +  jl'(z)  = 


^oh 


jvr 


60k  h  1  ~  Tt  I  7  .  „  (2),.  . 

o  J  -5 — J  (k  a)H  '  '  k  a 

L*  OO  O  O 


+  2 


n 

OO 

<  1 

-  +2  y 

l 

+  2  z 

m=  1 

00 

y 

cos 

L 

m=  1 

7  2 1  (k 
m  0  < 

1 7Tz/h 


o  o  m 


(1.  34) 

where  I"  is  the  component  of  the  current  in  phase  and  I1  is  the  component  of 
the  current  in  phase  quadrature  with  the  driving  voltage  and  the  admittance 


Y  =  (-y-)  is 

z=0 


Y  = 


60k  h  ]  .  TT  T  \tt  (2),,  : 

0  1  J  —  Jo(koa)Ho  (koa) 


+  2 


VJU 

2 


,  y  I  (k  ay  )K  (k  ) 
1  'm  o  o  'm  o  o  'm 


0  <  k  h 
o 


when 


<  7 J 


(1.35) 
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When  the  plate  separation  is  greater  than  a  half-wavelength 


From  1.. 35  it  is  seen  that  the  conductance  of  the  antenna  is  due  only  to  the 
primary  TEM  wave  corresponding  tom  =  0,  whereas  all  the  evanescent 
higher  modes  contribute  only  to  the  susceptance.  When  the  plate  separation 
is  increased  beyond  a  half-wavelength,  the  first  higher-order  modem  =  1 
starts  propagating  and  introduces  a  conductance  term,  as  can  be  seen  from 
1. 36.  Because  of  this  coupling  between  the  two  modes,  the  conductance  of 
the  antenna  is  now  larger,  signifying  increased  radiation.  As  the  plate 
separation  increases  further,  the  higher-order  modes  cause  a  periodic  increase 
in  the  antenna  conductance.. 

Closed-Form  Solutions  for  the  CurrentDistribution 
For  sufficiently  large  values  of  m 

7m  =  ^/(m  7T  Aoh)2  -  1  «  m  7T  /kQh 

and 

k  a.y  ~  m  7T  — . 
o  'm  —  h 
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Therefore,  for  the  higher-order  terms  the  I  and  Kq  functions  may  be  replaced 
by  their  asymptotic  values  given  by 

I  (x  )  e  /(27Tx  and  Kq  (x)  ~  (tt  /  2x)  e  X  for  large  x,  (1.  37) 

With  these  results  it  follows  that 


m  7T 


7  I  k  a>  K  k  a7  ~  ~  - 

'm  o  o  m  o  o  m  —  2k  a  —  ,  ,  , 

o  2  (k  a)  (k  h) 


-.(1.38) 


To  obtain  closed-form  expressions  for  the  current  the  series  in  1 ..  34 
can  be  truncated  at  M  h/ffa  and  .1  .  36  can  be  substituted  in  the  coefficients 
of  the  higher-order  terms. 


The  current  is  then  given  by 


j  5-  J  (k  a)H  (2)(k  a) 

J  2  o'  o  o  '  o 


+  2 


cos  m 


M 

— ,  y  u  I  (k  ay  ).K  (k  a-y  ) 
m=l  'm  o'  o  ‘m  o'  o  'm 


7T  z/h 


4(k  a)  .  (k  h) 
o  o 


OD 

l 


cos  m  n  z/h 


(1.  39) 


m=M 

The  series  in  the  last  term  of  the  right-hand  side  can  be  evaluated  by 
means  of  the  summation  formula 


UU 

I 

m=l 


cos  m 


W  z/h 


=  -  In 


2  sin  z/h 


With  this  formula  in  1.  39  the  current  distribution  is  given  by 


«■>  ■  TiFT 


_ 1 

j  -4U  J  (k  a)  .  H  (2)  (k  a) 
J  2  o  o  o  'o' 


(1.40) 
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2(k  a),  (k  h)  . 

'  o  '  o  '  1 


in- 1  m  o  o  m  o  o  m 


cos  mltz 


-  Aft  (kQa).  (kQh).lnJ  2  sin  (it  z/Zh) 


for  0  <  k  h  <  7T  . 
o 


(1.41) 


In  Fig.  1.  5  the  in-phase  component  I"  and  the  phase -quadrature 
component  I!  of  the  current  calculated  from  Eq.  1.  41  for  an  antenna  with 
a/X.  =  0.  01058  and  kQh  =  ir  /2  (h  =  X./4)  with  M  =  10  have  been  plotted.  The 
in-phase  component  I",  which  is  contributed  only  by  the  principal  radiating 
TEM  wave,  is  found  to  be  constant  over  the  entire  length  of  the  antenna;  it 
makes  a  significant  contribution  to  the  current  distribution.  The  phase- 
quadrature  component  I1  is  contributed  by  the  evanescent  higher-order  modes; 
it  involves  the  logarithmic  term,  and  hence  has  a  singularity  near  the  driving 
point  giving  rise  to  an  infinite  gap  susceptance;  this  aspect  will  be  discussed  in 
greater  detail  later.  The  logarithmic  term  actually  represents  the  current 
charging  the  adjacent  knife  edges  (with  infinite  capacitance)  of  the  idealized 
delta -function  generator.  I1  makes  only  a  minor  contribution  near  the  end, 
where  the  behavior  of  the  current  distribution  is  determined  mainly  by  I". 

f — - i - 

The  current  amplitude  I1  =J  I'  +  I"  therefore  tends  to  flatten  off  at 

Z  I  Z  Z 

the  end  of  the  antenna  as  can  be  seen  from  Fig.  1.  6  . 


The  current  distribution  for  antenna  lengths  in  the  range  it  <  kQh  <  2ir 


is  given  by 
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In  Fig.  1.  8  the  theoretical  values  (with  M  =  10)  of  I1  and  I"  for  an 
antenna  with  kQh  =  4.  7124  (h  =  3/4K),  a/\  =  0.  01058  have  been  plotted  . 
Now  the  contribution  to  the  in-phase  component  I"  comes  from  both  the 
TEM  wave  and  the  first-order  mode  (m  =1).  Since  the  first-order  mode 
has  a  cosinusoidal  dependence  on  length,  the  component  I"  no  longer  is 
a  constant  as  in  the  case  of  the  quarter-wavelength  antenna  (kQh  =  7T  /2). 

It  has  a  slightly  higher  amplitude  near  the  driving  point  than  at  the  end. 

The  I'  component  again  shows  a  logarithmic  rise  near  the  driving 
point;  it  has  a  minimum  at  k^z  =  1.  6  but  once  again  rises  to  a  maximum 
near  the  end. 
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Numerical  Accuracy  of  Calculations 

The  term  M  at  which  the  series  can  be  truncated  is  roughly  given 

by  h/tra  and  depends  on  the  thickness  and  height  of  the  antenna.  For  shorter 

and  thicker  antennas  fewer  terms  are  required.  In  practice,  however, 

the  number  of  terms  actually  needed  depends  entirely  on  the  desired  numerical 

tolerance,  For  the  case  of  the  exDerimental  antenna  with  k  a  =  0.  0664  and  with 

o 

k^h  ranging  from  0,  5  to  6.1a  summation  of  the  first  ten  terms  was  found  quite 

adequate  to  obtain  the  required  degree  of  accuracy  for  comparison  with 

experimental  results.  In  Fig.  1.4  l/lQ  (^0aTm)  •  Kq  (k^ay^)  has  been  plotted 

against  its  asymptotic  equivalent  2kQaym  for  an  antenna  with  h  =  X/4  or 

kQh  =  it /2  and  a/X  =  0.01058  and  also  for  a/X  =  0.03175.  The  term 

M  =  h/ffa  in  the  two  cases  is  approximately  8  and  3,  It  can  be  seen  from 

the  curves  that  the  error  in  the  truncation  process  is  indeed  very  small. 

Comparison  with  Experimental  Results 

In  Fig.  1.  6  and  Fig.  1.  7  the  theoretical  and  experimental  results 

r-2 — 

for  the  amplitude  of  the  current  |l|  =  y  I"  +  I'  and  the  phase  angle 

0  =  tan  ^  I'/l"  have  been  plotted.  As  can  be  seen  there  is  excellent 

agreement  between  the  two  except  near  the  driving  point  because  of  the 
logarithmic  singularity. 

In  Figs.  1.  8  and  1.  9  the  amplitude  and  phase  of  the  current  for  an 

antenna  with  kQh  =  - ^  “  ^X/4)  and  a/X  =  0.  010  58  have  been  calculated 

from  Eq.  1.  42  with  M  =  10,  and  compared  with  experimental  results. 

Agreement  is  very  good. 


0  4  8  12  16  20  24  28  32  36  40  44  48 
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FIG  1-4  CURVES  ILLUSTRATING  THE  VALIDITY  OF  THE  TRUNCATION  PROCESS. 
Io(koaym)  'Ko(koa  ym)  2  koaym  =  F2<m) 


MILLIAMPERES  PER  VOLT 


SERIES)  k0h=7772;  k0a  =  00664.  (a/X  =  0.0105); 
M  =  10. 


3.8 


FIG.  1-6  COMPARISON  OF  THEORETICAL  AND  EXPERI¬ 
MENTALLY  DETERMINED  CURRENT  DISTRIBU¬ 
TION  (AMPLITUDE)  k0hr7 r/2;  <VX*0;0105;  M*  10 


— » 

FIG.  1-7  COMPARISON  OF  THEORETICAL  AND  EXPERI¬ 
MENTALLY  DETERMINED  CURRENT  DISTRIBU¬ 
TION  (PHASE).  k0h  =  ir/2-,  a/x=0.0l05;  M=IO 


FIG.  1-8  THEORETICAL  CURRENT  DISTRIBUTION  (FOURIER 

SERIES  METHOD).  k0h=37r/2;  a /X  =  0.01058  ;  M  =  10. 


NORMALIZED  PHASE  ANGLE  (DEGREES)  - -  NORMALIZED  CURRENT  AMPLITUDE  1 1 1 


FIG.  1-9  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTALLY 
DETERMINED  CURRENT  DISTRIBUTION  (AMPLITUDE). 
k0h  =  37772  ;  a/X  *  0.0105  ;  M  » 10  . 


k„z  — - 


FIG.  1-10  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTALLY 

DETERMINED  CURRENT  DISTRIBUTION  (PHASE).  k0h  =  3ir/2; 
o/X  =  0.01058  ;  M  =  10. 
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Current  Distribution  and  Admittance  at  Resonance 

When  the  antenna  is  a  half-wavelength,  i.  e.  ,  kQh  =  7T ,  the  coefficient 

of  the  cosirz/h  in  Eq.  1.  42  shows  a  singularity  and  needs  close  inspection. 

(2)  2 

For  small  arguments  Jq  (x)  1  and  Hq  (x)  1  -  j  log  (  (3x  ) 

where  p  =  7  -  log  2  and  7  =  Euler's  constant  =  0.  5772.  With  these  identities 
used  in  the  coefficient  of  the  cosffz/h  term  in  Eq.  1.  36,  the  term  containing 
the  singularity  may  be  denoted  as  follows: 


When  kQh  =  it  (v  =  0),  P  =  cd.  This  gives  rise  to  an  infinite  current  on  the 
antenna.  This,  of  course,  is  physically  quite  meaningless  since  the  current 
on  the  antenna  has  always  to  be  a  finite  quantity.  A  similar  anomaly  has 
also  been  reported  by  Lewin  (3).  To  circumvent  this  difficulty  it  can  be 
assumed  that  the  parallel-plate  region  has  a  small  but  finite  attenuation 
constant  owing  to  losses  in  the  two  ground  planes,  so  that  k  is  now  a 
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complex  quantity  with  a  small  negative  imaginary  part  -  ja  ,  where  a  is 
the  attenuation  constant.  This  assumption  is  practically  quite  plausible 
since  the  image  planes  are  never  perfectly  conducting.  The  complex  wave 
number  can  be  represented  by  £  =  p-ja,  where  p  =  2ir/X. 

For  ph  =  7T 


If  the  a.^  term  is  neglected,  it  follows  that  v  =  -  ^  -  •  (1.  45) 

Also  £a  kQa.  Substitution  of  these  results  in  1.  43  gives  P  =  r/a 

(1.46) 


where  t 


(1.47) 
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Since  the  attenuation  constant  is  generally  a  very  small  quantity,  the 
cos  — ~  term  completely  predominates  over  all  other  terms  in  1,  42 
(except  for  the  logarithmic  term  close  to  the  driving  point)  and  so  the 
current  on  the  half-wave  resonant  antenna  becomes  a  pure  cosine 
I(z)  =  Pcos7T2j/h  (1.48) 

where  P  is  given  by  1.  46. 

This  phenomenon  has  been  confirmed  experimentally.  In  Fig.  1.11 
and  Fig.  1.  12  the  experimentally  observed  amplitude  and  phase  of  the 
current  on  the  antenna  have  been  compared  with  the  theoretical  curve.  There 
is  close  agreement  between  the  two.  The  admittance  of  the  half-wavelength 
antenna  is 

fY,.k  h=»  =  ?  -  ■  »• «) 

o 

This  is  a  very  large  quantity  which  is  limited  only  by  the  small  but  finite 
attenuation  constant  of  the  parallel-plate  region  to  the  propagating  waves. 

In  Fig.  1.  13  the  theoretically  calculated  values  for  the  conductance  and 
susceptance  of  the  antenna  in  the  vicinity  of  resonance  have  been  compared 
with  the  experimentally  measured  values.  In  Table  1-1  the  measured  values 
of  conductance  and  susceptance  of  the  antenna  near  resonance  have  been 
tabulated.  From  Fig.  1.13  and  also  Figs.  1.17  and  1.18  it  can  be  seen 
that  both  the  conductance  and  susceptance  show  a  sharp,  almost  discontinuous 
rise  near  k^h  =  7T  .  The  conductance  particularly,  as  can  be  seen  from 
Fig.  1.17,  shows  a  sharp  'spike'  when  the  antenna  is  nearly  a  half-wavelength 
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long.  The  conductance  jumps  from  2.  29  millimhos  at  k^h  =  3.  07  to  29.  64 

millimhos  at  k  h  =  3.  14  .  Since  the  conductance  is  very  sensitive  to  the 
o 

antenna  length  near  resonance,  the  experimental  measurements  had  to  be 
done  with  great  care. 

The  antenna  near  resonance,  therefore,  behaves  essentially  like  a 
short-circuited  transmission  line.  The  current  is  almost  a  pure  cosine 
and  its  amplitude  is  very  high.  An  examination  of  Figs..  1.  11  and  1*12 
shows  that  the  current  in  the  lower  and  upper  halves  of  the  antenna, while 
equal  in  magnitude,  are  180°  out  of  phase  and  so  their  individual  contributions 
to  the  field  at  distant  points  cancel  out.  This  implies  very  little  radiation 
loss  or  low  radiation  resistance  and  high  conductance. 

The  same  phenomenon  repeats  when  kQh  =  2 ir ,  only  this  time  the 
singularity  is  m  the  m=2  term  and  the  current  distribution  is  given  by 

I(z  )  oc  cos  2  ir  z/h  .  (1.  50) 

In  Fig.-  1,14  the  amplitude  of  the  experimentally  observed  curve  has 
been  compared  with  theory  to  confirm  the  above  result.  The  conductance 
again  shows  a  sharp  rise  as  can  be  seen  from  Fig.  1-17.  . 

A  more  rigorous  treatment  of  this  phenomenon  has  recently  been 
made  by  Tai  Wu  (6),  without  resorting  to  the  expedient  of  a  complex 
propagation  constant. 
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k  h 
o 

G[mhos  x  10^] 

3 

jB[mhos  x  10  ] 

2.  41 

1.13 

+j8.  02 

2.  63 

1.  00 

-  +jl0.  83 

2.  71 

1.  12 

+jll. 76 

2.  83 

1.  06 

+jl4,  72 

2.  89 

1.  16 

+jl7. 75 

2.  93 

1.  22 

+j20. 18 

2.99 

1.  28 

+j26.  65 

3.  02 

1.  16 

+j31 .  22 

3.  05 

1. 87 

+j 36 ,  92 

3.  07 

2.  29 

4j44. 14 

3.  11 

4.  18 

+j72. 77 

3.  14 

29.  64 

+jl89. 25 

3.  20 

18.  10 

-  j  79 . 02 

3.  30 

12.  68 

-  j  27 . 87 

3.  46 

6.  98 

-J9.81 

3.  68 

4.  20 

-j3.  62 

TABLE  1-1.  Measured  Values  of  the  Admittances  Near  Resonance 


(a  A  =  0.01058). 


FIG.  I-II  COMPARISON  OF  THEORETICAL  AND 

EXPERIMENTALLY  DETERMINED  CURRENT 
DISTRIBUTION  CURVES.  (AMPLITUDE). 
k,h  =  rr  ;  a/X  =  0.0105  ;  M  =  10. 


FIG.  1-12  COMPARISON  OF  THEORETICAL  AND  EXPERI- 
MENTALLY  DETERMINED  CURRENT  DISTRIBUTION 
(PHASE).  k„  h  =  TT  ;  o/X>  0.0105 


— h -G  4  I  Bl  x  ICr  mhos 


RELATIVE  CURRENT  AMPLITUDE 


FIG.  1-14  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTALLY  DETERMINED 
CURRENT  DISTRIBUTION  (AMPLITUDE);  k0h  =  27r;  a/X  =  0.0105. 
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Effect  of  the  Logarithmic  Singularity  in  the  Antenna  Current 

The  current  as  given  in  Eq.  1.  41  can  be  separated  into  two  parts 
I(z)  =  I  A(z)  +  I  (z)  ,  I.  (z)  represents  the  current  determined  by  the 

A  OO  A 

impedance  of  the  antenna  itself  and  is  given  by 


ZA  = 


JV_ 

6  Ok  h 


M 


j  ~  J  (k  a).  H 
J  2  o'  o  o 


+  2 


(kQa) 


II 


Ay  I  (k  ay  ).  K  (k  a  7  ) 
1  \ 'm  o  o  'm  o  o  'm 


2(k  a)(k  h) 
o  o 

7 T  m 


cos 


mf  z 


h 


(1.  51) 


On  the  other  hand, 

IQ0(Z)  -  -  -  ^  (kQa)  •  In  ^2  sin  (7Tz/2h)^  (1.  52) 

represents  the  current  determined  by  the  infinite  capacitance  of  the  knife 
edges  of  the  idealized  delta -function  generator;  it  has  a  singularity  at  the 
feed  point  giving  rise  to  an  infinite  input  admittance.  It  is  to  be  noted  that 

V 

only  the  phase -quadrature  component  of  the  current  is  infinite;  this  means 
that  only  the  input  susceptance  is  infinite  while  the  input  conductance  remains 
finite. 


The  amplitudes  of  |l|,  |l.  |  and  III  for  an  antenna  with  k  h  =  v  /Z 

A  CO  O 

and  a/\  =  0.  01058  have  been  plotted  in  Fig.  1. 16. 

Such  a  singularity  in  the  feed-point  current  has  been  investigated 
by  several  authors  (7,  8,  9,  10).  The  results  obtained  here  are  somewhat  at 
variance  with  the  corresponding  case  for  free-space  antenna.  The 
following  observations  may  be  made. 


SRI  5 


-28- 


1.  According  to  King  and  Wu((7)  and  Duncan  (8)  the  logarithmic  term  in  the 
case  of  a  free-space  antenna  is 


(z)=  -  ■ 


jV 


ooA  '  30tt 

For  the  parallel -plate  antenna 


(k  a)  .  In  k  z 
o  o 


IooP^Z'  =  "  lln  (koa)  ln  ^ 2  sm  ( nz/2h ) 


1 


(1.  53) 


(1.  52a) 


The  factor  2  missing  in  1.  52  arises  from  the  fact  that  a  monopole  is 
treated  instead  of  a  dipole  as  in  (7).  In  the  parallel -plate  case  the  logarithmic 
term  depends  on  the  height  of  the  antenna  (spacing  between  plates)  whereas  in 
the  free-space  case  it  is  independent  of  the  length  of  the  antenna.  This 
difference  is  probably  due  to  the  imaging  effect  of  the  two  plates  which  gives 
an  infinite  series  of  generators  separated  by  a  distance  h. 

2.  Though  such  a  singularity  in  the  feed-point  current  has  been  predicted 
theoretically  for  the  free-space  antenna,  no  such  phenomenon  has  been 
detected  by  either  the  Fourier  or  iterative  solutions  of  the  integral  equation 
for  the  antenna  current.  Duncan  and  Hinchey  (10)  have  carried  out  the 
Fourier  expansion  to  as  high  as  the  25th  order  without  noticing  any  untoward 
behavior. 

3,  Many  of  the  previous  workers  have  defined  a  'transition  region1  or 
'boundary  layer  region1  where  the  logarithmic  function  makes  a  significant 
contribution.  According  to  Wu  and  King  (7)  the  transition  zone  is  given  by 


z 

s 


=  \/Ze 


-l/Ka 


(1.  54) 
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For  the  antenna  used  in  the  experiments  k^a  =  0.0664  (a/X.  =  0.01058); 

_  7 

this  gives  k  z  =  2.  9  x  10 

6  os 

According  to  Duncan  (8)  the  corresponding  transition  zone  for  an 
antenna  with  kQa  =  0.  01  is  kQZg  =  0.  008  (zg/X  =  0.  0013).  He  has  also 
shown  that  zg  <  \/25  for  h/a  60.  In  Fig.  1.15,  the  amplitude  of  the 
current  near  the  driving  point  measured  experimentally  for  an  antenna  with 
k^h  =  7i  /Z  and  a/X.  =  0.  01058  has  been  compared  with  the  theoretical 
results  obtained  by  the  Fourier-s eries  method.  If  the  transition  zone  is 
defined  as  the  point  at  which  the  logarithmic  term  exceeds  the  nominal  value 
of  the  observed  current,  then  from  Fig.  1.15  the  transition  region  is  kQzg  = 
0.175  ,  which  is  much  higher  than  the  values  predicted  for  the  free-space 
antenna . 

4.  It  has  been  further  proposed  that  in  order  to  get  a  finite  value  for  the 
input  admittance  the  logarithmic  term  in  the  current  distribution  should 
be  subtracted  out.  Such  a  subtraction  process,  while  valid  for  thin 
structures,  may  be  questionable  for  thick  antennas  like  the  one  used  in  the 
present  experiment.  In  Fig,  1.16,  the  amplitude  of  the  currents 
|l|,  1 1^  |  and  j  1^  |  for  an  antenna  with  koh  =  7f/2  and  a/X.  =  0.010  58 
has  been  plotted.  It  can  be  seen  that  the  logarithmic  term  makes  a  marked 
contribution  over  a  significant  length  of  the  antenna.  In  fact,  the  dip  in  the 
current  distribution  at  k^z  =  0.2  seems  to  be  due  principally  to  the 
logarithmic  term.  However,  as  the  antenna  gets  thinner  the  logarithmic 
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term  plays  a  less  significant  role  because  of  the  factor  k^a  in  the  coefficient; 
in  fact,  for  an  infinitesimally  thin  antenna  it  has  been  shown  (see  ensuing 
discussion)  that  the  current  distribution  becomes  a  flattened  cosine  curve. 
Admittance  of  the  Antenna 

To  compare  the  theoretical  values  of  the  admittance  with  the 
experimental  results  we  have  to  subtract  out  the  infinite  gap  capacitance 
due  to  the  idealized  delta-function  generator.  Thus,  from  equation  1.  41 


However,  it  is  necessary  to  add  to  this  result  a  correction  term  to  account 
for  the  input  susceptance  of  the  actual  gap  or  feed  (9,  10).  The  correction 
term  denoted  here  by  j®q^pmay  be  obtained  empirically  by  comparing  the 
experimental  results  with  the  theoretical  values  for  B  given  in  equation 
1.  54.  The  final  theoretical  results  to  be  compared  with  experiment  are  then 

yf  =  Go  *  i(Bo  +  bgap>  ■=  °F  +  iBF  «•  55> 

Note  G  =  Gr. 

O  r 

In  Tables  1-2  (A)  and  1-2  (B)  the  gap  correction  ™  ^ 


o 


(1.  56) 


CURRENT  AMPLITUDE  — MILLIAMPERES  PER  VOLT  ^  - -  NORMALIZED  CURRENT  AMPLITUDE  (I) 


i.  1-15  COMPARISON  OF  THEORETICAL  AND 

EXPERIMENTALLY  DETERMINED  CURRENT 
DISTRIBUTION  NEAR  DRIVING  POINT, 
k.h  =  v/Z  ;  a/\  =  0.01058  ;  M  =  10 


FIG.  1-16  THEORETICAL  CURRENT  DISTRIBUTION  CURVE 
WITH  AND  WITHOUT  LOGARITHMIC  TERM. 
k,h  ■  TT/Z  ;  a/X- 0.0105;  M*IO. 
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BGAP  BEXP  ‘  Bo 

BEXP  =  Experimentally  Measured  Admittance 

B  =  Susceptance  After  Subtracting  Infinite  Gap  Susceptance  Obtained  by 
Fourier-Series  Method. 

=  Gap  Correction  Susceptance 
A) 


k  h 

0 

B^fmhos  xlO^] 

BEXP^mhoS  xl°3J 

BGApfmhos  xl°3l 

0. 7854 

-j 5.  49 

-jl.  10 

1. 0000 

-j3. 83 

+j0 .  20 

1. 2960 

-j2.  25 

+jl.  60 

1. 5708 

-j0.99 

+j3.  10 

1. 8450 

+j0.  21 

+j4.  50 

2. 0000 

+j  1 . 06 

+j  5.  30 

2. 3562 

+j3. 34 

+j7.  60 

Average  value  of  =  +j4.  10  for  0  <  kQh  <  ir  . 

B) 


k  h 

0 

B^fmhosxlO^] 

Bp^pirnnhos  x  10 

BGAP^mk°s  x  ^  ] 

3. 5000 

-j 12.  41 

-J7.40 

+ j  5 . 01 

3. 9270 

-J6.23 

-j0.  40 

+j5„  83 

4. 3000 

-J3..98 

+jl. 10 

+j5.  08 

4. 6200 

-  j-2.  55 

+j2.  80 

+j5.  35 

4. 7124 

-j2.  17 

+j3.  20 

+j  5.  37 

4. 9000 

-jl. 37 

+J4.  15 

+j5.  52 

5.  1000 

+j0 . 47 

+j5.  40 

+j4.  92 

5.  4978 

+j  2.  00 

+j5.  55 

+j5.  55 

Average  value  of  B^-^p  =  +j5.  41  Tl  <  kQh  <  Zv  . 

TABLE  1-2  (A)  and  (B).  Evaluation  of  Gap  Correction  Susceptance. 
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ha  s  been  tabulated.  *s  t'*ie  va-*-ue  susceptance  obtained  experimentally 

and  Bq  is  given  by  Eq.  I.  54. 

Because  of  their  high  values  and  steep  variation, the  values  of  B 

near  resonance  have  not  been  taken  into  account  while  calculating  B^^p  lest 

they  should  lead  to  erroneous  results. 

It  can  be  seen  that  in  the  range  0  <  kQh  <  7T ,  Is  nearly  constant 

_3 

and  has  an  average  value  of  4.  10  x  10  mhos.  In  the  range  7T  <  k^h  <  2tf  , 

_3 

B^-.^p  is  again  nearly  constant  but  average  value  is  5.  41  x  10  mhos.  This 
increase  in  B^p  correction  term  is  probably  due  to  the  fact  that  when  the 
plate  separation  is  greater  than  a  half-wavelength  ,  there  are  two  modes 
propagating  in  the  guide  which  might  change  the  field  distribution  at  the 
aperture  of  the  driving  coaxial  line. 

Comparison  with  Experimental  Results 

In  Figure  1-17  the  theoretical  values  of  the  conductance  obtained  by 
the  Fourier— series  method  have  been  compared  with  experimental  results. 


TABLE  1-3 


Gp  +  jBp  - 


B  + 
o 


BGAP 


Founer-Series  Admittance  of  Antenna 


where  Bq  =  Susceptance  after  Subtracting  out  Infinite  Gap  Susceptance 
Obtained  by  Fourier-Series  Method.  =  Empirical  Gap  Correction 

Sus  ceptance . 

B^p  =  +j4.  102  for  0  <  k^h  <  u 

B^p  =  +j5.  413  for  7T  <  k^h  <  2tt 
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where  D  is  an  infinite  positive  constant.  It  is  worth  while  noting  that,  as 
in  the  free-space  case  (7),  is  proportional  to  the  frequency  and 
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EB 

Gpfmhos  xlO^] 

i'B  [mhos  x  10^] 

J  o 

j'Bp[mhos  xlO^] 

0.  5000 

5.  00 

-  j8 .  8  3 

-j4. 73 

0. 7854 

3.  18 

-j5.  39 

-jl  •  29 

1.0000 

2.  50 

-j3.  83 

+j0.  26 

1. 2960 

1. 93 

-j2.  25 

+jl ■ 85 

1.  5708 

1.  59 

-j0.  99 

+j3. 10 

1. 8450 

1.  35 

+j0.  21 

+j4.  31 

2. 0000 

1. 25 

+jl.  06 

+j5. 16 

2. 3562 

1.06 

+j3. 34 

+J7.45 

2. 7000 

0.  92 

+j8.  27 

+j  1  2.  37 

2. 8250 

0.  88 

+jl  2.  01 

+ j  1 6 . 11 

3.1000 

0.  80 

+j8  2.  72 

+j86.  82 

3.  1300 

0.  79 

+j309. 47 

4- j  3 1  3.  57 

3.  1500 

96.  89 

-j334. 23 

-j 328 .  81 

3.  1700 

36.  29 

- j 1 10 . 84 

-jl 05.  43 

3. 2500 

13.  02 

-  j  34 .  10 

- j 28 .  69 

3, 3000 

9.  88 

-j-24.  64 

-  j 1 9 . 23 

3. 3570 

7.  94 

-  j  1 9 . 01 

-j 1 3.  59 

3.  5000 

5.  60 

-j  1  2.  41 

-j7. 00 

3. 9270 

3.  36 

-  j 6 . 23 

-j0. 81 

4. 3000 

2.  62 

-j3.  98 

+jl.  42 

4. 6200 

2.  25 

-j2.  55 

+J2.85 

4. 7124 

2.  16 

-j2.  17 

+j3.  24 

4. 9000 

2.  02 

-jl.  37 

+j4. 04 

5. 1000 

1. 88 

+j0. 47 

+j5.  88 

5. 4978 

1. 67 

+j2.  00 

+J7.41 

5.  8000 

1.  54 

+j7.  68 

+jl 3.  09 

6.  1500 

1. 42 

j - 

-  j  2  3 .  46 

-j  18.  05 

TABLE  1-3 


Theoretical  Admittances  of  Antenna  Determined  by  Fourier- 
Series  Method;  a/X.  =  0,  01058  ;  Truncation  term  M  =  10  . 


CONDUCTANCE  [mhos 


FIG.  1-17  COMPARISON  OF  THEORETICAL  AND  EXPERIMENTALLY 
MEASURED  CONDUCTANCE.  a/X  =  0.0105  ;  M  =  10 
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circumference  of  the  antenna  but  is  independent  of  the  antenna  length  h. 

Since  £  =  \/ \x  /£  and  k  =  u  \Jjj.  £ 

o  v  o'  o  o  v  o  o 

the  gap  capacitance 

Y 

Coo  =  =  8  CQ  a  D  =  8  So  a  Lt  [  y  +  In  M]  .  (1.  62) 

J  M-^oo 

It  has  already  been  noted  that  a  finite  gap  correction  term 

must  be  added  to  the  theoretical  susceptance  values  in  order  to  bring  them 

into  agreement  with  experimental  results.  For  0  <  k^h  <  IT  the 

correction  term  evaluated  empirically  from  experimental  results  is 

4. 10  x  10  ^  mhos . 

~  -jL-  .  (k^a)  .  [  y  +  In  m  ]  ~  4, 10  x  10  ^  mhos . 

The  corresponding  value  of  m  =  10.  74  ~  10. 

ThuSj  the  gap  correction  to  be  added  to  Yq  in  Eq.  1,  54  is 


B 


GAP 


j/60  k  h. 
o 


4(kQa)  (kQh) 


10 


1 

m 


=  4.1x10' 


mhos . 


m=  1 


This  is  exactly  the  last  term  on  the  right-hand  side  of  Eq.  1.  54  which  was 
subtracted  out.  This  means  that  in  Eq.  1.  3  5  we  need  add  only  the  first  10 
terms,  which  represents  the  contribution  to  the  susceptance  from  the  first  10 
higher-order  modes,  to  get  results  for  the  susceptance  which  agree  with 
the  experimentally  measured  values. 
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Current  Distribution  on  an  Infinitesimally  Thin  Antenna 

The  general  expression  for  the  current  distribution  on  the  antenna  is 
given  by  1.  34' .  Consider  now  the  term  in  this  equation  which  involves  the 
series  summation  and  denote  it  by  S  . 

oo 

S  =  2  T  -  cos  muz/h -  (L  63) 

,  7  .  I  (k  a7  ).  K  (k  ay  ) 

m=l  ‘m  o  o  'm  o  o  m 

For  small  arguments,  K^(x)  =  log  (px),  and  Iq(x)  ~  1,  where  log  (3  = 

y  -  log  2  ~  -  0. 1159, 

/  /  m7T  N\  ,  m7T 

j{~ iTh— )  -  1  ~  -ETT-  • 

Hence  , 

=(-£v! — )  •  ” 

and 

iog  (PkQa7m)  =  log  |  +  log  m  .  (1.  64) 

We  note  that  as  a-*  0,  log  ( j3a7T  /h)  -»  oo. 

This  means  that  for  an  infinitesimally  thin  antenna,  with  vanishing  radius, 
the  log  ((3a7T  /h)  term  completely  predominates  over  the  log  m  term. 
Therefore , 


iog  (pkoa7m)  ~  log  (pa7T  /h) 


(1.  65) 
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(1.69) 


SRI  5 


-36- 


This  is  of  the  form  l(:z)  =  A  +  B  cos  kQ  (h  -  z),  and  the  current  distribution 
is  a  slightly  flattened  cosine  curve. 

Current  Distribution  with  a  Filamentary  Kerntel 

In  the  solution  of  the  integral  equation  for  the  antenna  current  using 
iterative  methods,  it  is  customary  to  replace  the  correct  kernel  distance 
of  the  form  1,  10  by  an  approximate  'filamentary'  kernel  distance  (1,  8) 
given  by 


Rm  =V('Z  1  2'-2mh)2  +  a2  .  (1.70) 

It  is  assumed  that  h  >  >  a. 

While  such  an  approximation  may  be  valid  for  the  free-space  antenna 
problems,  it  leads  to  divergence  difficulties  in  the  current  in  the  case  of 
wave-guide  probe  investigations.  Thus, if  1.  70  xs  substituted  in  1.  5  and 
the  manipulations  carried  out  the  current  distribution  on  the  antenna  comes 
out  as 


jV 


60k  h 
o 


i  \  Ho'2'<V> 


+  2 


7T  z/h 


m=l 


y  .  K  (k  ay  ) 
'm  o  o  'm 


•(1.  71) 


For  higher  values  of  m  ,  y  ~  mu  /k  h  and 
°  ’  m  —  '  o 

,  U  a 

k  a-v  ~  - = —  m  . 

o  'm  —  h 


The  asymptotic  expansion  for  ^Q(x)  Is"" 


(IT  e  ”X 


therefore.  l/K  (k  ay  )  .  y 
>  '  o  o  'm  'm 


2 


behaves  like 
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2/<ko*>  ■  <koh>3  • 

Thus, the  series  for  the  current  distribution  diverges.  Such  a  divergence 
difficulty  has  been  noticed  by  Lewm  (4)  while  investigating  the  problem 
of  a  probe  in  a  waveguide. 

Distribution  of  Charge  on  the  Antenna 

The  charge  distribution  may  be  ob  lined  directly  from  the  equation 
of  continuity 

dl 

-■az— '  +  J  wclz  =  0  (L72) 

where  I  =1"  +  jl'  is  given  by  Eq.  1.  34  .  Thus  from  1.  72  the 

z  z  z 

distribution  of  charge 


a+(7 T  a/h).  m 
3/2 


81" 


81' 


q(z)  =  j(q"z  +  jq'z) 


j 

u 


+  j 


so  that  for  0  <  k  h  <  j 
o 

,  81" 


8  z 
8l',z 


0.  Since  I"(z)  is  a  constant  for  k^h  <7 T  .  (1.  73) 


8  z 


30u 


<koh> 


Y"  ms 

Z  tJ.u 


msinm 


IT  zi/h 


m=  1 


o'(V\n)Ko(W 


The  charge  distribution,  being  sinusoidal,  is  maximum  near  the  center  and 

vanishes  at  the  ends.  When  k  h  >  77,  q"  4  0  because  of  the  contribution 

o  ’  ^  z  ' 

from  the  first  higher-order  mode  m  ~  1. 
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3.  VARIATIONAL  METHOD 

It  has  been  shown  in  Section  2  that  while  the  Fourier-series  method 
gives  accurate  solutions  for  the  current  distribution  and  conductance  of 
the  antenna,  the  values  for  the  susceptance  obtained  from  it  are  ambiguous 
because  of  the  feed-point  singularity  in  the  current.  The  variational 
technique,  however,  like  the  one  employed  by  J.  Storer  (11)  in  the  case  of 
a  free-space  antenna,  ignores  the  singularity  of  the  driving  point  and  provides 
an  approximate  but  continuous  solution  for  the  current  distribution  over  the 
entire  length  of  the  antenna. 

The  stationary  expression  for  impedance  can  be  formulated  in  the 
following  manner: 

h 

v0  6  (ss)  =  j  ?Q  J°  K(z,  z ')  I  (z1)  dz1  . 
o 

The  driving  voltage  Vq  -  Iq  where  is  the  input  impedance  and  I 
input  current.  Therefore, 

h 

ZoIo5(z)  =  j?o  f  K(z,  z')  I  (*:')  dz'  .  (1.75) 

o 

If  both  sides  are  multiplied  by  I(z)  and  integrated  between  0  and  h  and  then 
2 

divided  through  by  I  (0)  the  result  is 

h 

J '  K(z,  z')  I(z )  I  (z1)  dz  dz1  (1.  76) 

o 
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where  K(z,  z1)  is  the  kernel  given  by  1.  29  and  1.  30.  Since  the  kernel  is 
symmetric,  i.e.,  K(z,  z')  =  K(z',  z),  Eq,  1.  76  can  be  shown  to  have  all  the 

variational  characteristics,  and  Zq  is  stationary  with  respect  to  small 
changes  in  I  . 

According  to  1.  34  and  1.  69,  a  suitable  trial  function  for  the  current  is 
IT  (z)  =  A  4  B  cos  k  (h  -  z)  .  (1.77) 

With  this  type  of  current  distribution  the  impedance  is  given  by 


dz  dz 1 


(1.80) 
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h  h 

Vgg  =  J*  J  K(  z,  z')  cos  kQ(h  -  z).  cos  kQ(h  -  z')  dz  dz'  .  (1.  81) 

o  o 

The  two  complex  coefficients  for  the  trial  current  A  and  B  have  to  be 
determined  in  order  to  calculate  the  impedance  from  1.  78..  It  is 
assumed  that  a  unit  voltage  is  driving  the  antenna;  the  current  at  z  =  0 
is  determined  by  the  following  condition! 


V  =  l  =  IT(0)  .  Zy  . 


I  (0)  =  A  +  B  cos  k  h  =  l/Z, 


(1.82) 


F rom  (1.  82) 


B  cos  k  h 
o 


V 


(1.83) 


ZV  j?o  (Z)2  YAA 


B  cos  k  h 
o 


+  2"ab-  b 


-  B  cos  k  h 


JV 


vBB  •  B  )>  •  (1.84) 

Zy  is  an  analytical  function  of  the  parameter  B.  Thus,  the  next  step  is  to 
8Z 

equate  —  to  zero  and  solve  for  the  optimum  value  of  B  which  makes 
2iy  stationary. 

We  obtain 


B  = 


VAA  COS  koh  -  VAB. 


(1.85) 


JV 
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where  A  =  v..  cos  k  h  -  2d._  cosk  h  +  Vr,„  •  (1.86) 

AA  o  Ax5  o  no 


Using  (1.  85)  and  (1.  83)  we  get 


A  = 


r  v .  „  cos  k  h 
BB  AB  o_ 


ZV  A 


(1.87) 


If  '1.87'  and  1.85'  are  substituted  in  1.78  the  result  is 

2 


Z  =  j? 
V  J5,o 


_VAA'  VBB  "  VAB  J 


(1.88) 


Evaluation  of  the  'v1  integrals 

The  various  v  integrals  occurring  in  the  above  equations  can 
now  be  evaluated;  the  following  simple  integrals  will  prove  useful  in 
the  derivation: 


J 1  cos  kQ(h  -  z)  dz  =  sin  kQh  /k 


and 


j '  sin  kQ  (h  -  z)  cos  ™  ^  Z  dz 


sin  k  h 
o 


1  -  (m7T  Aph)^ 


(1.89) 


From  Eq.  1.  79  we  have 
h  h 


_  j  Jr^ 

vAA  ~  f  f  K(z,z')  dz  dz'  =  - f—  Jo(koa)  Ho(2)(koa).  (1.90) 


o  o 
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From  (1.  .80; 


h  h 


AB 


z,z')  cos  kQ(h  -  z)  dz  dz1 


o  o 


sin  k  h  ... 

j  — Jo  <V>  ■  Ho  <koa| 


(1.91) 


0  <  k  h  <  2u  . 
o 


From  1.8-1! 

h  h 


BB 


=  f  K(Z,zl)  cos  kQ(h-z)  cos  kQ  (h-z1)  dz  dz' 


o  o 


-  1 


27rk  h 
o 


j  \  J  (k  a)  .  H  (2)  (k  a). 
J  2  o  '  o  o  '  o 


sin  k  h 
o 


+  2 


oo 


I 


I  (k  a  y  ' )  .  K  (k  ay  )  , 

o'  o  o  o  'm  .  2 


sin  k  h  V  for  0  <  k  h  <  v  . 


m=l  rm 

For  higher-order  terms 


I  (k  ay  )  .  K  (k  ay  )  ~  — == - 

o  o  'm  o  o  'm  —  2k  ay 

o  m 


k  h 
o 


2(kQa).  m77 


Therefore 
oo 


J 

m-M 


I  (k  ay  ).K  (k  ay  )  (k  h)' 
o  o  'm  o  o  'm  o 


2  (k  a) .  7T 
o 


M 


m  •  ,  m 

m-i  m=l 


where  M  is  the  truncation  term  in  the  series. 
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<koh)3 

2.  (kQa)  .  7T  3 
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+  2 


M 

I 

m=2 


I  (k  a  7  ).K  jk  av  ) 

o  o  'm  o  o  'm 


(koh> 

2i? .  k  a 
o 


+  1.  202 


(koh>3 

V  3. (kQa) 


(1.93) 


The  numerical  evaluation  of  the  admittance  Yy  =  l/Zy  =  Gy  +  jBy 
and  the  current  distribution  parameters  A  and  B  has  been  carried  out 
for  k^a  =  0.0664  (a/X.  =  0.01058)  with  M  =  10.  The  numerical  values 
of  A  and  B  are  given  in  Table  1-4  and  they  are  plotted  in  Fig.  1-20. 

From  Figures  1-17,  1-18,  and  1-19,  it  is  seen  that  the  theoretical 
admittance  values  are  in  good  agreement  with  experimental  results  until 
kQh  =  3.  7  but  deteriorate  rapidly  thereafter.  The  theoretical  admittance 
values  are  tabulated  in  Table  1-5.  Thus,  the  trial  function  for  the  current 
of  the  form. A  +  B  cos  kQ  (h  -  z)  is  valid  only  in  the  range  0  <  k^h  <  3.  7 
but  a  modification  will  probably  have  to  be  made  for  longer  lengths,  to 
account  for  the  first-order  mode  which  propagates  when  kQh  >  n  .  In 
Table  1-6,  the  theoretical  value  for  the  susceptance  B^^  obtained  by 
the  variational  method  are  compared  with  the  susceptance  values  B0 
obtained  by  the  Fourier-series  method.  The  difference  -  BQ  is 

almost  constant  up  to  kQh  =  2,  8250  and  the  average  value  of  By0  = 

+j3.  63  ■  this  compares  favorably  with  the  gap  correction  term  Bq^p  = 
+j4.  10  that  has  to  be  applied  to  BQ  obtained  by  the  Fourier- series  method. 
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The  average  difference  |B^J  =  I  ®EXP  ”  *s  only  |  0 .  3  5  J  millimhos 

indicating  good  agreement.  In  Table  1-7  susceptance  obtained  by 

the  variational  method  has  been  compared  with  experimental  results. 

The  current  distribution  for  an  antenna  with  k  h  =  7T  /2  and  k  a  =  0.  066 

o  o 

has  been  calculated  from  1.  77  and  the  results  are  compared  with  the 
experimental  curve  in  Fig.  1.  21.  It  can  be  seen  that  the  general 
agreement  is  not  as  good  as  it  was  in  the  case  of  the  Fourier-series 
solution  Fig.  1-6,  In  Fig.  1-15  the  current  near  the  driving  point  has  been 
compared  with  the  experimental  curve  and  the  Fourier-series  solution. 

Though  the  qualitative  agreement  with  experiment  is  rather  poor,  the 
variational  method  does  not  show  any  abrupt  'logarithmic'  rise  near  the 
feeding  point  as  in  the  case  of  the  Fourier-series  solution.  It  is  interesting 
to  note  that  it  is  possible  to  obtain  the  far -field  radiation  pattern  of  the  antenna 
in  terms  of  A  and  B. 

Impedance  of  a  Very  Short  Antenna 


When  the  plate  separation  is  small  we  may  assume  the  current  on 


the  antenna  to  be  a  constant.  If  we  also  assume  a  thin  antenna, then 


jk  ? 

J  oso 


h  h 


iFTT1  •  1  S  f  j  "T"  Jo(koa)-  Ho<2>(koa)  dz  dz' 

^  o  o 


UU  XJL  a* 

-2  V  I  (k  ).  K  (k  a7  )  y  2  P  P 
/_j  o  o  m  o  o  m  'm  J  J 


cos  mir  z  cos  m7T  z 1 

h  ‘  h 


dz  dz 


(1-94) 


SRI  5 


-46- 


Z  =  30 77  (k  h)  J  (k  a)  .  H  (2)  (k  a) 
o  o  o  o  o 


For  k  a  «  1  J  (k  a)  ~  1 
o  ,  o  o 


H  (2)  (k  a)  =  1  -  j  .  —  {-  0. 1159  +  log  k  a 

O  O  J  7T  \  6  O 


Hence,  the  impedance  of  the  antenna 

Z  =  3077  (kQh)  +  j  60  kQh  .  log  (l/(3kQa) 


Cl.  95) 


The  above  result  is  the  familiar  formula  for  the  impedance  per 
length  h  of  an  infinitely  long  wire  carrying  a  uniform  current. 


k  h 
o 

A  x  10  3 

B  x  103 

0. 5000 

0, 496 

+j3. 900 

-0. 368 

- j4 1 .  221 

0. 7854 

3.  366 

+jl  2.  293 

-0. 239 

-  j  20 .  105 

1. 0000 

17. 334 

+jl0.  336 

-17.  643 

- j 1 7 .  650 

1.  2960 

1.948 

+j4,  292 

-0. 031 

-jlO. 446 

1.  5708 

1. 650 

+j 2 .  776 

-0. 011 

-j  8. 873 

1. 8450 

1.  361 

+jl.  826 

+0.012 

- j  8 . 195 

2. 0000 

1.  254 

+jl.  457 

-0. 003 

- j  8 . 161 

2.  3562 

1. 064 

+j0 .  813 

-0. 004 

-j  9.  085 

2. 7000 

0.  907 

+j0.  315 

+0.  528 

-jll, 919 

2.  8250 

0.  753 

+j0.  553 

+7. 094 

-  j  1 4.  850 

3.  3570 

0.  315 

-jO.  344 

-6. 584 

+j  1  5.  703 

3. 5000 

0.  262 

-jO.  335 

-4.  525 

+j9  .  510 

3.  9270 

0.  207 

-jO.  495 

-2.  392 

.+j 3  .  625 

4.  3000 

0.  316 

-jO. 814 

-1 .  253 

+jl  .  185 

4. 6200 

0.  854 

-jO. 018 

+1.  504 

-jO  .  440 

4. 7124 

1.  154 

-jO.  993 

+  2.  931 

-jO  .  173 

4. 9000 

1.  374 

-jO.  360 

+4.  306 

+j2 . 784 

5.  1000 

1.  234 

+j0. 030 

+4. 091 

+j5  .  310 

5. 4978 

0.  702 

+j0.  207 

+0.993 

+j5  ..686 

5.  8000 

0.  594 

+j0.  252 

+2. 335 

+j  1  3.  200 

TABLE  -  1-4  . 


Current  Distribution  Parameters 
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Yy  =  Gy  +  jBy  =  Variational  Admittance 


k  h 
o 

Gy[mhos 

x  10  3] 

3 

Byfmhos  x  10  ] 

0. 5000 

4. 

71 

-5.  89 

0. 7854 

3. 

18 

-1.76 

1. 2960 

1. 

94 

+  1. 45 

1. 5708 

1. 

60 

+  2.  77 

1. 8450 

1. 

35 

+4.  04 

2. 0000 

1. 

26 

+4.  85 

2.  3562 

1. 

06 

+7.  23 

2. 7000 

0. 

92 

'  +11.44 

2.  8250 

0. 

89 

+  15.  42 

3.  3570 

6. 

88 

-15.  69 

3. 5000 

4. 

52 

-9.  25 

3.  9270 

1. 

89 

-3.  06 

4.  3000 

0. 

81 

-1.  27 

4. 6200 

0. 

72 

-1.03 

4. 7124 

1. 

15 

-0.99 

4.  9000 

2. 

17 

+0.  16 

5.  1000 

2. 

78 

+  2.  08 

5.4978 

2. 

05 

+  5.  60 

5.  8000 

2. 

66 

+11. 94 

TABLE  1-5.  Theoretical  Admittance  of  Antenna  Determined  by 

Variational  Method;  a/X.  =  0.01058  ;  Truncation  TermM=10 . 
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B,,  =  B...  _  -  B 

Vo  VAR  o 

B^r  =  Susceptance  Obtained  by  Variational  Method 

B  =  Susceptance  after  Subtracting  Infinite  Gap  Capacitance  Obtained 
o'  " 

by  Fourier-Series  Method. 


k  h 
o 

B.r  =  Blr .  _  -  B, 

V  o  VAR  o 

0. 7854 

+ j 3 .  63 

1.  2960 

+j3.  70 

1.  5708 

+j3.  76 

1. 8450 

+J3.  79 

2. 0000 

+j3.  79 

2. 3562 

+ j 3 .  89 

2. 7000 

+j3.  17 

2.  8250 

+ j  3 .  31 

Average  value  of  Byo  =  +  j 3 .  63 

TABLE  1-6.  Comparison  of  Susceptance  By^^ Obtained  by  the  Variational 
Method  with  the  Susceptance  Bq  after  Subtracting  Infinite 
Gap  Susceptance  Obtained  by  Fourier-Series  Method. 
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ba  bexp  '  bvar 

BVAR  =  ®uscePtance  Obtained  by  Variational  Method 
BEXP  =  ^easureb  Susceptance 


k  h 
o 

3 

BVAR[mhosxl°  ] 

BExp[mhos  xl°3] 

bexp~bvar 

=Bjmhos  x  10  3] 

0=.  78.54 

— j  1 . 76 

-jl.  10 

+j0 . 66 

1. 2960 

+ j  1.45 

+jl.  60 

+j0.  15 

1.  5708 

+j2.  77 

+j3. 10 

+j0. 33 

1. 8450 

+j4.  04 

+j4.  50 

+j0.  46 

2. 0000 

+j4.  85 

+j  5.  30 

+j0 .  45 

2.  3562 

+j7. 23 

+j7.  60 

+j0.  37 

2. 7000 

+jll. 44 

+jll.  50 

+j0.  06 

2.  8250 

+jl  5.  42 

+j  1  5.  00 

-j0.42 

3.  3570 

4jl 5.  69 

-jl  5.  50 

+j0. 19 

Average  value  of  |  |  =  0.  35  millxmhos 


TABLE  1-7.  Comparison  of  Susceptance  Obtained  by  the  Variational 

Method  with  the  Experimentally  Measured  Susceptance  Bgxp 
(a/X.  =  0.01058;  M  =  10). 


FIG.  1-21  COMPARISON  OF  THEORETICAL  (VARIATIONAL) 
AND  EXPERIMENTAL  CURRENT  DISTRIBUTION; 
k0h  =  7T/2;  V0.0I05;  M  =  10 
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4,  EXPERIMENTAL  EQUIPMENT  AND  MEASURING  TECHNIQUES 

This  section  is  devoted  to  the  experimental  phase  of  this  research. 

The  experimental  arrangement  can  be  divided  into  two  major  groups: 

1.  The  parallel-plate  region  with  the  associated  absorbing  wedges  for 
reducing  edge  reflections;  2,  The  coaxial  measuring  line  and  ancillary 
equipment  for  measuring  the  input  admittance  and  current  distribution 
of  the  antenna. 

Choice  of  Frequency 

The  experimental  design  of  the  equipment  was  complicated  by  the 
problem  of  selecting  a  suitable  frequency  which  was  a  compromise  between 
two  contradictory  factors  - 

1.  It  was  not  practically  convenient  to  have  ground  planes 
larger  than  6  ft.  square  because  of  the  size  limitations 
of  the  laboratory  room;  hence,  the  frequency  had  to  be 
high  enough  so  that  the  error  due  to  the  finite  size  of  the 
ground  planes  was  not  very  significant.  Besides,  inordinately 
large  and  heavy  ground  planes  were  undesirable  since  it 

would  then  be  very  difficult  to  vary  the  plate  separation 

i 

without  resorting  to  complicated  mechanical  devices. 

Moreover,  as  the  frequency  decreases  it  becomes  increasingly 
difficult  to  find  suitable  absorbing  material  to  terminate  the 
rims  of  the  parallel-plate  region  to  minimize  edge  reflections. 


/ 
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2.  The  frequency  at  the  same  time  had  to  be  sufficiently  low 
in  order  not  to  limit  the  accuracy  of  the  measuring 
equipment,  because  of  the  mechanical  tolerances  involved 
and  to  minimize  deleterious  effects  like  probe  loading  and 
large  end-corrections.  Furthermore,  the  plate  spacing 
had  to  be  varied  by  small  fractions  of  a  wavelength  in  order 
to  study  the  behavior  of  the  input  admittance  of  the  antenna 
as  a  function  of  its  length;  this  could  be  accomplished  more 
easily  at  longer  wavelengths. 

It  was  proposed  to  carry  out  the  measurements  at  a  frequency  of 
1000  Mcs.  (L,  Band)  which  seemed  to  be  a  reasonable  compromise  between 
these  two  extremes  and  had  the  additional  advantage  that  convenient 
oscillators  were  readily  available. 

Design  and  Construction  of  the  Parallel- Plate  Region 

The  error  due  to  the  plates  being  finite  instead  of  infinite  as 
theoretically  required  is  difficult  to  estimate  as  no  adequate  experimental 
data  are  available.  However,  from  a  study  of  the  results  obtained  for 
antenna  impedance  measurements  in  free  space  it  was  decided  that  the 
minimum  size  of  the  plates  should  be  six  wavelengths  to  obtain  reasonably 
accurate  results.  Meier  and  Summers  (13)  have  studied  the  effect  of  screens 
smaller  than  6X  on  the  measured  input  admittance  of  antennas  in  free  space. 
For  a  square  ground  plane  of  side  6k  they  have  reported  a  maximum  error 
of  10  for  a  thin  antenna  of  length  h/X  =  0.  224  and  an  error  of  +  25  ohms 
for  a  length  of  h/X  =  0.448.  The  percent  error  in  |  Z  |  could  thus  be 
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almost  3  °/o  ;  the  errors  are  found  to  be  slowly  damped  oscillatory 

functions  of  the  dimensions  of  the  ground  screen.  The  corresponding  errors 

in  the  case  of  a  parallel -plate  region  are  likely  to  be  of  a  lower  order  of 

magnitude  since  the  edge  reflections  can  be  virtually  eliminated  by  using 

absorbing  wedges..  To  reduce  the  error  further  the  antenna  was  mounted 

off-center  by  a  quarter -wavelength.  Andrews(14)  reports  that  for  an  antenna 

-2 

placed  off-center  the  error  varies  as  d  where  d  represents  the  side 
dimension  of  the  image  plane. 

In  accordance  with  the  above  requirements  the  over-all  dimensions 

of  the  parallel-plate  device  were  made  6  feet  x  6  feet  (or  approximately 

6  wavelengths  square  at  1000  Me  ).  Actually  the  bottom  plate  measured 

1" 

6  feet  x  12  feet  (6\  x  12X)  and  was  fabricated  from  —  Dural  sheets 

mounted  on  a  sturdy  wooden  framework  to  prevent  sagging.  To  enable  the 

measurement  of  the  standing-wave  ratio  inside  the  parallel-plate  region  the 

1" 

lower  plate  was  divided  lengthwise  by  a  sliding  panel  made  from  — 

Dural  sheet  6  ~  feet  long  x  1  inch  wide.  The  panel  carried  a  small  electric 
probe  near  its  center  and  could  be  driven  by  means  of  a  geared-down  motor 
and  a  cable  sliding  over  a  pulley.  The  position  of  the  probe  could  be 
determined  to  0.1  mm  by  means  of  a  steel  scale  and  a  vernier  mounted 

I 

on  the  lower  plate.  Fig.  1.  22  shows  a  picture  of  the  sliding  panel  j 

probe  and  driving  mechanism.  The  top  plate  was  deliberately  made 

very  light  so  that  the  plate  spacing  could  be  varied  easily  by  merely 

1  " 

raising  or  lowering  the  upper  plate.  It  was  fabricated  from  a 
Dural  sheet  mounted  on  a  light  plywood  frame.  The  spacing  between 
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two  plates  was  maintained  by  wooden  posts  placed  at  the  four  corners; 

these  seem  to  have  little  effect  on  the  field  distribution  within,  since 

the  field  strengths  at  the  edges  of  the  plates  are  relatively  low.  The  plate 

separation  could  be  varied  by  small  distances  with  the  help  of  thin 

plastic  slabs  placed  on  the  supporting  wooden  posts.  During  the  course  of 

1  " 

the  experiment  the  plate  spacing  (antenna  height)  was  varied  from  — 

(0.  04X.)  to  12"  (1.  016X).  While  examining  the  behavior  of  the  antenna  near 
resonance  the  plate  spacing  had  sometimes  to  be  varied  by  intervals  as  small 
as  0.  3  cm  (0.  01X).  The  plates  were  aligned  by  means  of  metal  pins  which 

fit  through  holes  at  the  sides  of  the  two  plates  so  as  to  index  the  upper  - 

plate  relative  to  the  lower;  after  adjusting  the  plate  spacing  by  means  of  the 
wooden  corner  posts,  the  metal  pins  were  removed  to  avoid  possible 
perturbing  effects  during  the  course  of  the  measurements.  The  joints  and 
screw  slots  on  the  surface  of  the  plates  were  covered  with  strips  of  .  001" 
aluminum  foil  with  silver  paint  as  adhesive.  The  two  plates  were  carefully 
levelled  to  within  0,  010  inch  with  the  help  of  a  theodolite  and  any  unevenness 
noticed  was  corrected  by  placing  suitable  shiims  beneath  the  plates.  Row  (15) 
has  investigated  the  effects  of  small  variations  in  plate  spacing  and  it  is 
found  that  in  practice,  gradual  departure  from  the  ideal  has  very  little  effect 
on  the  performance,  though  sudden  discontinuities  may  cause  partial 
reflections.  Bracewell  (16)  has  found  that  the  effect  of  these  discontinuities 
is  to  introduce  a  shunt  capacitance  (  a  )  in  the  parallel-plate  line  whose 
magnitude  is  given  by 
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C 


d 


I 


.  2 

sin  n7T  a 

”2  3 


a  n 


where  a  is  the  ratio  of  the  spacings  (see  Fig.  t-23  )and  Cq  (.  085 
picofarads/meter)  is  the  permittivity  of  free  space.  If  the  discontinuity 
is  not  large  (a;  >,0.9)  the  discontinuity  capacitance  is  only  a  minute 
fraction  of  a  picofarad  and  does  not  cause  noticeable  perturbation.  The 
other  references  which  were  helpful  in  designing  the  parallel-plate  region 
were  (17 -20)  •  ' 

Absorbers 

The  wave  propagating  inside  the  parallel-plate  region  suffers  a 
partial  reflection  at  the  edges  because  of  the  discontinuity  in  impedance. 
The  theoretical  value  of  the  reflection  coefficient  as  a  function  of  plate 
spacing  has  been  calculated  by  Marcuvitz  (21)  and  is  plotted  in  Fig.  1-24  . 
It  can  be  seen  that  strong  reflections  occur  for  close  plate  spacing  which 
can  seriously  affect  the  measured  value  of  the  admittance  of  the  antenna. 

In  order  to  eliminate  the  reflection  it  was  necessary  to  terminate  the  rims 
of  the  parallel  plate  with  absorbing  wedges.  The  wedges  had  to  be 
sufficiently  thin  (approximately  1  cm)  so  that  the  plate  separation  could 
be  varied  in  small  intervals  in  order  to  study  the  behavior  of  the  antenna 
admittance  as  a  function  of  its  length.  Commercially  available  L-band 
absorbers,  besides  being  prohibitively  expensive^  were  also  far  too  large 
to  be  suitable  for  the  purpose  of  this  experiment.  It  was  therefore  decided 
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to  make  the  absorbers  in  the  laboratory  itself  from  linearly  tapered 
polyfoam  wedges  coated  with  an  energy  dissipating  substance  like  carbon 
black  or  graphite.  The  two  determining  factors  (22)  in  the  design  of  these 
'gradual  transition'  type  absorbers  are  (1)  the  length  of  the  transition, 

(2)  the  optimum  layer  thickness  of  the  absorbing  substance.  The  lower 
frequency  limit  is  given  by  the  length  of  the  wedges.  Halford  (23)  has  shown 
that  the  reflection  coefficient  of  such  wedges  varies  periodically  with  their 
length  as  shown  in  Fig.  1,  25  .  There  appears  to  be  no  clear-cut  upper 
limiting  frequency  though  in  the  case  of  acoustic  absorbers  (24),  it  is 
found  that  the  reflection  rises  at  higher  frequencies  when  separation 
between  the  individual  wedges  becomes  too  large.  Because  of  its  low 
dielectric  constant  and  the  ease  with  which  it  could  be  cut  and  tapered 
the  wedges  were  made  from  —-inch  thick  polyfoam  sheet  and  were 
tapered  for  1  wavelength  (11  3/4  inches).  The  separation  between  the 
individual  wedges  was  approximately  5  cm  (  X)  and  the  total  length  of 
each  wafer  was  4  feet.  Thin  plywood  strips  provided  with  wooden  screws 
were  glued  to  the  base  of  these  wedges.  The  wedges  were  then  sprayed 
with  3-E-14D  Black  No.  1  --  a  colloidal  dispersion  of  carbon  black  in  liquid 
Neoprene,  a  product  manufactured  by  the  General  Latex  and  Chemical  Co.  , 
Cambridge,  Mass.  The  whole  wafer  thus  assembled  (see  Fig.  1.  26) 
could  be  inserted  inside  the  parallel-plate  region  and  screwed  onto  slotted 
wooden  panels  placed  at  the  edges  of  the  parallel  plates.  The  plate 


FIG-  1-25  VARIATION  OF  REFLECTION  COEFFICIENT  WITH 
LENGTH  FOR  LINEAR  TAPER.  (FROM  REF  23) 


FIG.  1-26  ABSORBING  WEDGES 


1 _ i _ I _ i _ I 

/  INCH 

*  > 


FIG.  1-27  SHIELDED  LOOP  CURRENT  PROBE 


28  PARALLEL- PLATE  REGION  WITH  TOP-PLATE  REMOVED 


FIG.  1-29  SIDE  VIEW  OF  PARALLEL- PLATE  REGION 


NORMALIZED  AMPLITUDE 


FIG.  1-30  FIELD  DISTRIBUTION  INSIDE 
PARALLEL- PLATE  REGION 
WITH  l  WITHOUT  ABSORBERS 
k0h  =0.5319 


SRI  5 


-53- 


separation  could  be  easily  varied  by  removing  or  adding  each  wafer  in  turn. 
Fig.  1.  28  shows  a  picture  of  the  parallel-plate  region  with  the  top  plate 
removed;  the  absorbers  mounted  on  wooden  panels,  the  antenna  and  the 
sliding  panel  with  probe  for  measuring  the  standing-wave  ratio  in  the  guide 
can  all  be  seen  in  the  picture.  Fig.  1.  29  gives  a  side  view  of  the  same 
arrangement  with  the  top  plate  in  place. 

The  favorable  layer  thickness  of  the  carbon  black  mixture  was 
determined  by  giving  the  wedges  various  coatings  and  determining  the 
standing-wave  ratio  in  each  case.  Fig.  1.  30  gives  the  normalized  field 
distribution  inside  the  parallel-plate  region  with  and  without  the  absorbing 
wedges  for  a  plate  separation  of  h/X.  =  0.  084  .  The  standing-wave  ratio 
without  absorbers  was  1.  72  at  a  distance  of  about  25  cm  from  the  antenna. 
After  the  insertion  of  the  absorbing  wedges  with  about  4  coatings  of-the  carbon 
black  mixture  the  SWR  was  brought  down  to  about  1. 12.  No  appreciable 
reduction  in  SWR  was  noticed  by  increasing  the  number  of  coatings  on  the 
wedges.  The  normal  component  of  the  electric  field  in  the  parallel-plate 
guide  is  given  by  Ez  =  cH^^  (kQr)  where  kQ  =  2tt/\.  This  has  also  been 
plotted  in  Fig.  1.  30  for  the  sake  of  comparison.  It  can  be  seen  from 
Fig.  1.  24  that  the  edge  reflections  become  less  pronounced  as  the  plate 
separation  increases.  This  fact  was  also  confirmed  experimentally.  For 
large  spacings  the  SWR  fell  to  less  than  1.  07  .  However,  when  the  plate 
separation  was  a  half-wavelength  a  sudden  increase  in  SWR  was  noticed. 
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The  reflection  coefficient  calculated  from  Marcuvitz's  admittance  values 
does  not  account  for  this  phenomenon.  After  inserting  the  absorbers, 
however,  the  SWR  was  once  again  brought  down  to  a.  reasonable  level. 
Measuring  Techniques 

The  experimental  techniques  for  measuring  the  impedance  and 

current  distribution  of  the  antenna  have  become  fairly  standardized,  thanks 

to  the  efforts  of  the  previous  research  workers  in  this  laboratory. 

The  coaxial  measuring  line  was  quite  similar  to  the  one  devised  by 

T.  Morita.  (25);  a  schematic  sketch  of  the  line  is  given  m  Fig.  1.  31.  The 

measuring  line  was  installed  vertically  beneath  the  parallel-plate  region. 

In  order  to  reduce  the  end  effects  and  the  possible  effect  that  the  aperture 

of  the  coaxial  line  might  have  on  the  multiple  series  of  images,  the  ratio 

of  the  diameters  of  the  outer  to  the  inner  conductor  was  chosen  to  be  2.  25  . 

1  " 

The  center  conductor  of  the  coaxial  line  was  — diameter  brass  tubing; 

it  could  be  extended  from  0  to  32  cm  above  the  lower  ground  screen  by 

1" 

means  of  a  rack  and  pinion  arrangement .  The  inner  conductor  had  a  ^  slot 

along  its  entire  length,  through  which  the  shielded  probe  could  travel. 

3  " 

The  probe  was  mounted  on  a  — jg —  brass  tube  which  slid  inside  the  center 
conductor  and  was  positioned  by  another  rack  and  pinion  arrangement. 

Both  the  probe  and  antenna  could  be  positioned  to  0. 1  mm  by  means  of  a  steel 
scale  and  vernier.  The  signal  from  the  probe  to  the  detector  was  carried 
by  a  microdot  cable.  The  coaxial  line  was  driven  by  means  of  a  T-shaped 
stub,  connected  to  the  main  line  at  a  distance  of  X/4  from  the  short- 
circuited  end.  The  probe  was  small  shielded  loop,  approximately  2.  5  mm 


FIG.  1-31  SKETCH  OF  MEASURING  LINE  AND  PARALLEL- PLATE 


REGION 


RELATIVE  CURRENT  AMPLITUDE 


FIG  1-32  SHIELDED  LOOP  PROSE  SYMMETRY  TEST  ON  SHORT 
CIRCUITED  LINE. 


FIG  1-33.  DISTRIBUTION  CURVE  IN  SHORT  CIRCUITED 
MEASURING  LINE.  (  TEST  FOR  PROBE 
LOADING) 
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(0.008X.)  in  diameter;  it  was  made  from  .021"  O.  D.  x  .003"  wall 
precision  coaxitube.  (See  Fig.  1.  27).  Since  the  frequency  of  1000  Me 
used  in  these  experiments  is  considerably  higher  than  frequencies  hitherto 
used  in  this  laboratory  for  antenna  measurements  of  this  nature,  a 
comprehensive  series  of  tests  was  made  to  determine  the  accuracy  of  the 
measuring  equipment.  Particular  attention  was  focused  on  the  preparation 
of  the  probe  which  usually  causes  the  greatest  errors.  Care  was  taken  to 
construct  the  probe  symmetrically,  since  in  the  vicinity  of  a  current 
minimum,  slight  asymmetries  can  introduce  sizeable  errors.  That  the 
final  probe  is  well-balanced  is  shown  by  the  symmetry  of  the  distribution 
curve  near  a  current  minimum,  (Fig.  1.  32).  In  order  to  determine  the 
possible  perturbing  effect  of  the  probe,  the  coaxial  line  was  short-circuited 
and  the  experimentally  measured  current  distribution  curve  inside  the  line 
was  compared  with  an  ideal  sinusoidal  distribution.  As  can  be  seen  from 
Fig.  1.  33  the  probe  does  cause  slight  perturbation  at  the  peak  of  the  curve, 
but  the  effect  does  not  appear  to  be  significant  enough  to  cause  appreciable 
errors.  The  input  impedance  of  a  free-space  monopole  as  measured  with 
the  present  experimental  set-up  compared  favorably  with  the  values 
obtained  by  Hartig  (27).  The  constants  of  the  measuring  line  are  listed  below: 
Inner  radius  of  outer  conductor  b  =  9/32  inch. 

Outer  radius  of  antenna  a  =  l/8  inch. 

b/a  =  2.25 
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Characteristic  admittance  of  line  =  Y  =  20.  57  mhos. 

c 

Characteristic  impedance  =  48.  6  ohms. 

Frequency  of  measurement  =  1000  Me. 

Wavelength  in  air  =  30  cm. 

Wavelength  in  line  =  28.  9  cm  ,  . 

Capacitive  end-correction  =  B^,  =  u  C^,  =  -0.  584  mi'llimhos . 

Apparatus 

A  superheterodyne  detecting  system  was  used  in  the  experiment  . 

The  signal  generator  used  in  the  experiment  was  a  General  Radio  Type 
1218-A  Unit  Oscillator  having  a  frequency  range  of  900-2000  Me  .  and  a 
maximum  power  output  of  200  MW.  Higher  harmonics  generated  by  the 
oscillator  were  suppressed  by  a  GR  Type  874-F1000  low-pass  filter  with 
a  cut-off  frequency  of  1000  Me.  (+  10  °/o  ).  To  avoid  frequency  pulling 
the  oscillator  was  isolated  from  the  rest  of  the  circuit  by  a  20  db 
attenuating  pad.  The  detecting  system  consists  of  a  General  Radio  Type 
1216  A-I.  F.  Amplifier  used  in  conjunction  with  a  GR  Type  874-MR  mixer 
rectifier  with  a  type  1-N21B  crystal  diode.  The  local  oscillator  was  another 
GR  Type  1218-A  Unit  Oscillator.  The  I..  F.  amplifier  is  provided  with  a 
70  db  precision  step  attenuator  with  an  accuracy  of  +  0.3db  +  l°/o  ; 
good  linearity  is  obtained  over  the  whole  range  of  the  instrument.  The 
amplifier  was  calibrated  by  using  the  sinusoidal  voltage  distribution  in  a 
short-circuited  coaxial  line  as  a  reference.  The  calibration  curve  is  shown 
in  Fig .  1.  34. 


FIG.  1-35  BLOCK  DIAGRAM  OF  EXPERIMENTAL  EQUIPMENT 
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A  schematic  sketch  of  the  experimental  layout  is  shown  in  Fig.  1.  35. 
Figure  1.  36  gives  a  pictorial  view  of  the  apparatus. 

Current  Distribution 

The  current  distribution  on  the  antenna  can  be  completely  described 
by  measuring  the  relative  amplitude  and  phase.  The  amplitude  can  be 
measured  quite  easily  by  connecting  the  probe  output  directly  to  the 
detecting  system  through  a  tuning  network. 

In  order  to  measure  the  relative  phase,  the  signal  from  the  probe 
is  compared  with  a  reference  signal  obtained  from  a  non-resonant  slotted 
line  terminated  in  its  characteristic  impedance.  The  phase  line  used  m 
the  experiment  was  a  Hewlett-Packard  Model  805A  slab  line  terminated 
by  a  GR  Type  874  WM  50  ohm  termination;  the  standing-wave  ratio  on  the 
line  was  less  than  1.  0  5.  The  probe  and  the  reference  signals  are  mixed 
in  by  a  coaxial  hybrid  junction  (Narda  Model  No.  3032).  The  difference  arm 
of  the  hybrid  was  connected  to  the  detecting  system.  By  adjusting  the 
phase  and  amplitude  of  the  reference  signal  from  the  phase  line  a 
maximum  indication  could  be  obtained,  signifying  that  the  two  signals 
are  180°  out  of  phase.  After  establishing  such  a  reference  point,  the  phase 
at  other  points  along  the  antenna  could  be  determined  with  respect  to  it 
by  simply  noting  the  distance  between  the  minimum  settings  of  the  probe 
in  the  phase  line. 
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Admittance  Measurement 

Two  techniques  were  employed  to  measure  the  input  admittance  of 
the  antenna,  At  low  VSWR  the  ratio  of  maximum  to  minimum  current  was 
measured  directly  and  the  terminal  function  p  of  the  antenna  load  was 
calculated  from  the  formula 


p  =  coth  1  S 
K  sa 

where  S  =  current  standing-wave  ratio.  At  higher  VSWR  the  terminal 
function  was  determined  by  the  distribution  curve  (double  minimum)  method. 
If  A  u  is  the  distance  between  half-power  points,  the  apparent  terminal 
function  is  given  by 


sa 


where  p  = 


where  \g  is  the  wavelength  in  the  coaxial  line. 

The  apparent  phase  function  can  be  found  from  the  equation 


$sa  =  "f-+  P  <W1  -  u2> 

where  is  the  position  of  the  current  minimum  with  a  short-circuit 
terminating  the  line  and  is  the  current  minimum  with  the  antenna 
load.  The  input  admittance  may  then  be  calculated  from  the  formulas. 


G _  sin  h  2  p _ 

Y  cos  h  2n  +  cos  2  $ 

c  r 

_  sin  2  { 

c  cosh  2 p+  cos  2  $ 
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and  Ya  =  G  +  jB  ;  =  charge  admittance  of  line  . 

End  Correction 

The  ideal  or  theoretical  antenna  admittance  Yq  may  be  obtained 

from  the  corresponding  measured  apparent  admittance  Y^  by  noting  that 

Y  =  Y  +  jUC 
a  o  J- 

where  C,p  is  a  lumped  terminal  zone  capacitance  that  takes  account  of  end 
and  coupling  effects  near  the  end  of  the  measuring  line.  C^,  has  been 
evaluated  approximately  (26)  and  is  given  by 

CT  =  -  f  (b/a)  cQb 

where  c^  =  Zn  SQ/ln  (b/a)  is  the  capacitance  per  unit  length  of  the  line,  a 
is  the  antenna  radius  and  b  is  the  inner  radius  of  the  outer  conductor  of 
the  coaxial  line. 

The  function  f(  —  )  =  -  C^/c^b  has  been  given  in  (26).  For 

b/a  =  2.25,  -  CT/cob  =  0.18937.  The  corresponding  end-correction  is 
-juC^,  =  -  0.  5839x10  mhos  .  The  end-correction  in  (26)  is  for  a /\  - 

-3 

3.  9x10  .  The  correction  to  be  applied  to  the  experimental  results  is  likely 

to  be  slightly  larger  since  a/k  =  10.  05x10  However,  this  difference  in 
correction  term  is  almost  negligible,  particularly  since  the  capacitance  end- 
correction  is  very  small  and  can  be  virtually  ignored  except  in  the  vicinity 


of  resonance  and  antiresonance. 
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APPENDIX  -  I 


Equation  1.  28  may  be  rigorously  derived  in  the  following  manner. 
Consider  the  integral 
7T 

^  (j2kQa  sin  Q)  dO  (1.  28) 


Q 


F rom  the  identity 


(j2kQ  a  sin  0)  =  -  j  7r / 2  [2k^a  s^n 


it  follows  that 
7T 


Q=-j7r/2  J  Hq^^[c  sin0]de  =  -  j  /2 


7 T 

I 


7/ 


J  (c  sin  G)  dO  -  j  \  N  (csin  0)  dO 


where  c  =  2k  a  , 
o 


(1.96) 


Using  the  formula 

n 

J ^  (z)  =  ^  Jq  (2z  sin  /)  cos  2  n  /  dj/  for  n  =  0,l,  2.  .  . 


we  have 


77 

Jq  (c  sin  0)  d0  =  7T  (c/2) 


(1.97) 


Consider  now  the  integral 
7T 


I 


N  (c  sin  9 )  d  0 
o 


o 
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The  integral  representation  for  the  N  function  is  given  by  5[M  and  O,  p.  27] 


N  (x)  =  - 


,  ,  1  00 

2  (  -y  x)  p 

Terror;  -  \ 


COS  X.  t 


r(f-v).r(T 


(t2-l)v  +  T 


dt 


for  -l/2  <  Re  v  <  -i-  ;  x  real  and  positive. 


When  v  =  0  we  have 


Since 


00 

-  2/*r  J 

J1 


77  00 


^  Nq  (c  sinO)  dO  =  -  2/77  ^  ^ 


cos  (c  sin  0)  t 


\ft‘ 


dt  dO 


o  1  \/ 1  -  1 

77 


d.t 


1 


cos  (c  t  sin  0)d  0 


t  -  1  o 


cos  (c  t  sin  0)  dO  =  77  Jq  (ct)  , 


(1.98) 


I  No 


(c  sin  0)  dO  =  -  2 


00 

I 


J  (ct) 


dt  . 


1  \/t~  -  1 


Making  the  substitution  t  =  \fp  ^  +  1,  we  have 


J  (ct) 
o 


dt 


UL> 

-f 


Jo  «= 


Jpz  +  1 


s/p  2  +  1  ) 

dp 


(1.99) 


(1. 100) 


Using  another  identity  given  by  (5) 
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f  j„ ,x\/‘r*  i) 

Jv/2(x/21  Nv/2  <*/2>  =  -2/"  7 , 


we  have  for  v  =  0 


dt 


/ 


t  +  i 


(Re  v  >-l,  x  real  and  positive) 


UU 

I 


j  (x/777) 


/.  2  7 


dt  =  -77  /2  J  (  x/2)  N  (x/2) 
'  o  o 


o  v  t  +  1 

Using  1.100  and  1 .  101  in  1 . 99 
77 


^  (c  sin  8)  d6  =77  Jq  (c/2)  Nq  (c/2) 


Substituting  1.  97  ,  L  102  in  1.  96  we  obtain 
77 

Q  =  ^  Kq  (j2kQa  sin  Q)  =  -  j  77  2 /2  Jy  (ka)  Hq(2)  (k(. 


(1.  101) 


(1.  102) 


(1.  28) 
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Closed  form  expressions  for  the 
current  distribution  and  impedance 
of  an  antenna  in  a  parallel-plate 
region  have  been  obtained  by  an 
integral  equation  technique.  Inter¬ 
esting  features  connected  with  the 
antenna  are  (1)  the  logarithmic 
singularity  in  the  current  near  the 
driving  point  because  of  the  ideal¬ 
ized  delta -function  generator  (2) 
behavior  at  resonance.  The  theoret¬ 
ical  results  have  been  verified  ex¬ 
perimentally  and  a  description  of  the 
apparatus  is  included  in  the  report. 
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